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This research concerns linear ordinary differential equations depending 
in such a way on a parameter p that the "limit" differential equation 
obtained by letting p tend to * in the differential equation is of lower 

\ r 

order than the original one. f t C 

v • 

Adopting a term customary in physics we used the name boundary layer 
problem for the question: What happens to the solution of a boundary value 
problem of such a differential equation, if the parameter tends to 00 in this 
solution? 

We gave a general answer to this question for the differential equation 
■^N(y) +M(y) » 0, where N(y) and M(y) are linear differential 
expressions of order n and m, respectively (n > m), and for non- 
homogeneous boundary conditions which consist in prescribing the values of 
derivatives (but not of linear combinations of such derivatives) at the end¬ 
points. The question whether the solution of such a boundary value problem 
converges to a solution of the limiting differential equations, as p ♦ 
and what boundary conditions are satisfied by the limit function could be 
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decided by an easily applicable rule. This rule showed, among other things 
that the solution converges only, if the prescribed n boundary conditions 
are not too unevenly distributed between the two endpoints. 

If the order ra of the limiting differential equation is only one less 
than the order n of the original differential equation, then the above 
mentioned rule could be extended to more general types of boundary conditions 
and also to non-homogeneous differential equations. 

Since the most important boundary layer problems in the applications are 
concerned with systems of differential equations, we gave a simple example for 
the mathematical treatment of a boundary layer problem for a linear system of 
two simultaneous differential equations. 

The validity of the general rule proved in this research was seen to be 
restricted by the assumption that the coefficient of the term of highest order 
of differentiation in M(y) has no zeros in the interval of integration. In a 
special example we showed that interesting results can be obtained, if this 
assumption is dropped. 

The theory of the asymptotic expansion of the solutions of linear 
differential equations involving a parameter, developed by G. D. Birkhoff, 
Noaillon, Tamarkin, Trjitzinsky and others proved an important and powerful 
tool in these investigations. 
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COMMENT 


In May 1980, the Mathematics Research Center organized a successful 
advanced Seminar on Singular Perturbations and Asymptotics in honor of the 
retirement of a colleague, Wolfgang R. Wasow. His fundamental research is 
responsible for many other rapid developments in this field since 1940, and 
continues to play a vital role in modern theory and current applications. 
Wasow's Ph.D. dissertation (N.Y.U., 1941), a small part of which exists in 
print (On the asymptotic solution of boundary value problems for ordinary 
differential equations containing a parameter, J. of Mathematics and Physics 
12 (1944), 173-183), represents the starting point of this Important flourish 
of modern applicable research. 

Following suggestions of several participants MRC is printing his 1941 
thesis in its entirety as a TSR in order to make this valuable work more 
widely available. Readers will note that the name "singular perturbations" 
(which was only coined several years later by K. 0. Friedrichs or w. Wasow or 
possibly jointly, but neither is now able to recall the details) does not 
appear anywhere explicitly! 
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ON BOUNDARY LAYER PROBLEMS IN THE THEORY OF 
ORDINARY DIFFERENTIAL EQUATIONS 

Wolfgang R. Wasow 
Introduction 

Many problems In applied mathematics lead to questions of the following type: 

Given Is a differential equation Involving a parameter 
p. This parameter occurs In such a way that the 
’limiting* differential equation, i.e. the differential 
equation obtained by letting p tend to Infinity In 
the differential equation , is of lower order than the 
original one. What happens then to the solution of a 
boundary value problem of the original differential 
equation, if p tends to infinity In that solution? 

It Is by no means obvious - and not even always true, as we shall see - that the 
solution of such a boundary value problem tends to a solution of the limiting differential 
equation, as p tends to Infinity, But even when this is the case the question arises 
what are the boundary conditions satisfied by the limiting function. As a solution of a 
differential equation of lower order than the original one it cannot, in general, be 
expected to satisfy all the boundary contions prescribed in the original problem. 

In those cases in which the solution of the original problem converges - as p tends 
to infinity - to a solution of the limiting differential equation which does no longer 
satisfy all the originally prescribed boundary conditions, the solution f the original 
problem shows a peculiar behavior for very large values of the parameter p. Some of the 
derivatives of the solution will assume very large values in a narrow region near the 
boundary. As p tends to infinity, these derivatives will tend to infinity at a certain 


The preparation of the report is sponsored by the United States Army under Contract No 
DAAG29-80-C-0041. 
















part of the boundary. In the most important applications of phenomena of this type the 
first derivative of the solution - and, of course, all the higher derivatives - diverge at 
parts of the boundary, as p tends to infinity. 

In the physical interpretations this means the occurrence of "Boundary layers" in 
which the quantity to be investigated increases or decreases very rapidly with the distance 
from the boundary, when some physical constant is large. We shall use the name Boundary 
Layer Problems , in a more general sense, for all related mathematical problems. 

The most famous problem of this type is the relationship between the theories of 
viscous and ideal liquids. An interesting boundary layer problem for a system of two non¬ 
linear ordinary differential equations has been investigated recently by X. Friedrichs and 
J. J. Stoker in a paper on the buckling of elastic plates, [9], 

The majority of the applications lead to non-linear partial differential equations 
which are so complicated that a complete mathematical treatment has not yet been attempted. 
But even the boundary layer problem for linear ordinary differential equations, a problem 
interesting from the mathematical as well as from the physical point of view, has as yet 
been hardly investigated. The only paper known to the author of this investigation, on 
this problem for ordinary differential equations is the article [7] by Erich Rothe, in 
which the problem is solved for a very special linear differential equation of the second 
order with constant coefficients. 

In Chapter I of the present paper we discuss the boundary layer problem for linear 
homogeneous differential equations depending linearly on the parameter, and for non- 
homogeneous boundary conditions. The result of this part can be expressed by a simple and 
easily applicable rule which determines immediately, for a given problem of this type, 
whether the solution converges and what boundary conditions are satisfied by the limiting 
function. 

In Chapter II we investigate more thoroughly the case where the order of the limiting 
differential equation is lower by one than the order of the original differential equation. 
In this case the statement of the general rule of the Main Theorem in Chapter I can be 
formulated so as to include more general boundary conditions than those assumed in Chapter 










I. The problem is then solved - at leact for a drop of one in the order of the 
differential equation - for the non-hoe jeneous equation. Finally, the usual method of 
treating boundary layer problems, consisting of a transformation of the independent 
variable and an appropriate modification of the boundary conditions, is shown to be 
justified in this case. The method is sometimes used in more complicated problems without 
mathematical justification. 

In Chapter III some special examples of other types of boundary layer problems are 
discussed. 

The methods used in this paper are based on the theory of the asymptotic solution of 
ordinary differential equations involving a parameter. This theory has been developed by 
G. D. Birkhoff [1], Noaillon [2], Tamarkin [3], [4], Trjitzinsky [6], and others. In the 
Appendix we give a short outline of the results of this theory as far as they are used in 
this investigation. 

I am deeply indebted to the Professors R. Courant and K. 0. Friedrichs 'hose help and 
encouragement played a major part in the preparation of this thesis. The original 
suggestion for this investigation came from Prof. Friedrichs, and his active interest in 
the progress of the work has been of the utmost value. 
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Chapter I 


THE BOUNDARY LAYER PROBLEM FOR THE DIFFERENTIAL 

EQUATION - N(y) + M(y) - 0 
P 

$1. Statement of the Problem 
1. He consider the ordinary linear differential equation 

L(y,p) « 0 , 

where the linear differential expression L(y,p) is of the form 

L(y,p> = - N(y) + M(y> 

P 


with 


( 101 ) 


( 102 ) 


n (n-v) 

N(y) = l a (x) y (x) (103) 

vO V 

m (m-y) 

M(y) = l b (X) y (x) . (104) 

P-0 U 

x is a real variable and p a positive parameter. He assume that the coefficients 

a (x) and b (x) admit at least n bounded derivatives in the interval 
\) u 

a « * < B . 

If the order n of the differential expression N(y) is greater than the order m 
of the differential expression M(y), then the differential equation (101) gives rise to a 
boundary layer problem for the "limiting" differential equation, i.e. the differential 
equation obtained by letting p tend to infinity in the original differential equation 
L(y,p) - 0. For, this limiting differential equation is 

M(y) - 0 , (105) 

and this differential equation is of lower order than (101), if 

n > m (106) 

He shall also assume that 

m > 0 . (106a) 


- 4 - 










Most of our results reaaln valid for a ■ 0. But at some points the Inclusion of the 
case m - 0 would make the statement of the result rather Involved. It seemed therefore 
preferable to exclude this case from the Main Theorem. 

Together with the differential equation (101) we prescribe n boundary conditions for 



the function 

y(x). The boundary conditions 

considered in 

this chapter 

are of the form 



L i <y> “ 4 i ' 

(i - 1,2. 

n) 

(107) 


with constant 

and with 






( <X i‘ 




1 


y <6) 

for i “ 1,2, 

r • « » /I 


4 


V y) * { , , 



(108) 







* 


l y fa) 

for i “ r+1. 

r e • » ,11 , 



where x - a and x » 8 are the left and right endpoints, respectively, of the interval 
under consideration.* 

We assume that the boundary conditions are arranged in such a way that 

X. > X, > ••• > 1 
1 2 r 


Vi > Va > — > T n • 

This arrangement is the opposite of the customary one, but it is more practical in our 

case. Ml the numbers 1^ and are, of course, assumed to be less than n. r is the 

number of boundary conditions prescribed at the right endpoint. The number of boundary 

conditions at the left end point is then n - r. 

One or both of the numbers X and t may be zero, which means that the value of 

r n J 

the function itself is prescribed at one or both endpoints. But our theory applies also to 
cases in which only derivatives of the function are prescribed at the endpoints. The 
boundary conditions (108) contain as a special case the initial value problem. We have 
only to set r “ 0, or r “ n. 


SI 


* d K 

we shall use throughout this paper the notation y* **(x) for — y , 

dx* 














We make further the assumption 


Sq(x) 9 < 0, for ail x in o < x <, 6 , (109) 

which makes it possible for us to set 

a 0 (x) - 1 . (110) 

without loss of generality. 

A very essential condition for the validity of the theory that follows is that we must 
have also 

b Q (x) f 0, for all x in a < x < 6 • (111) 

It is easy to see that a theory of boundary layer problems which does not assume (111) mist 
be expected to be much more complicated. For b Q (x) is the coefficient of the first term 
ii the limiting differential equation (105). Hence, if b Q (x) has zeros in a < x < B, 
these zeros will, in general, be singular points of the solutions of the limiting 
differential equation. 

To these assumptions we will have to add two more conditions of a rather essential 
nature. Since these conditions can be more easily formulated in connection with our Main 
Theorem we postpone their statement for a few pages. 

In general, the differential equation (101) will have a unique solution U(x,p) 
satisfying the boundary conditions (107). U(x,p) depends on the value of the parameter p. 
We will be able to give a general rule which allows us to decide, for a given problem, 
whether 

u(x) - lim U(x,p) (112) 

p+<» 

exists, and which are the boundary conditions satisfied by u(x). We shall see also that 
u(x>, when it exists, is a solution of the limiting differential equation M(y) “ 0. 

The behavior of U(x,p), as p tends to infinity, will be seen to depend, in 
general, on three things only: 

(a) On the number n - m, i.e. the difference between the orders of the 
original and the limiting differential equation. 
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(b) On r, i.e. on the way in which the n boundary conditions are 
divided between the two end points* 

(c) On the sign of the coefficient b Q (x)* 

2. There are a great many different possible cases for our boundary layer problem* In 
9ome cases U(x, p) converges, as p ♦ in some cases it diverges, and there are boom 
special occurrences that are not covered by the Main Theorem* This accounts for the fact 
that the Main Theorem, although very simple to apply, is somewhat lengthy to formulate* He 
precede its general formulation by a few examples, in order to give, without proof, an idea 
of the variety of possible occurrences* In the convergent cases the boundary conditions 
satisfied by the limit function u(x) are obtained by canceling n - m of the given 
boundary conditions, usually taken among those involving higher orders of differentiation* 
Example 1 . 

L(x,p) z — y"' + x 3 y + 2 x y" » 0 
P 

with the boundary conditions 

y" ( a) - l 3 y"(“) “ ^ 

y’lfl) * t 2 * 

If, e.g., a ■ 1, 8 “ 2, then b Q (x) >0 in a < x < 8 and the solution U(x,p) of the 
problem tends to the solution of the differential equation 

M(y) s x 3 y" ♦ 2 x y - 0 

satisfying the boundary conditions 

y" ( 8) - t, 
y’<8) - * 2 , 

which are obtained by canceling the boundary condition given at x “ a* If a ” -2, 

8 = -1, then b Q (x) < 0 and U(x,p) tends to the solution of M(y) • 0 with the 
boundary conditions 

-7- 


irri'- "* TT ~‘“ 1 " • * ^ 










I 


1 

* 


y-(a) - ^ 


yMB) - i. 


obtained by canceling the first boundary conditions at x « 0. 

But if a “ -I, B “ 1« the condition (111) is no longer satisfied and our Main 
Theorem does not apply• 


Example 2. 


L(x, p) = (y* 4 * + cos x»y* 3 * ) xy" + xy - 0 


with the boundary conditions 


y"<a) - i. 


y'(o) - i. 


If a “ -2, B “ -1. then b-(x) < 0 and 0(x, p) converges to the solution of the 


differential equation 


satisfying the boundary conditions 


y* + y ■ 0 


y* ( o) “ tj 
y (o) - 


which are obtained by canceling the first boundary condition at each endpoint. If a > 0, 
B > 0 then b Q (x) > 0 in a < x < 8, and U(x, p) tends in general, to the solution of 
y* + y » 0 with the boundary conditiras 


because in this case the general theory requires the canceling of the two boundary 
conditions involving the highest order of differentiation. 


a“ 2* , B“ , 

then we have an exceptional case. Because then there is no solution of y" ♦ y ■ 0 
satisfying the boundary conditions y( cO - i^, y( B) “ unless ^ » 0. Again, our 

Main Theorem does not cover these special values of a and B. 


Example 3. 


My, p) = 1 y <4) - xy' - y - 0 


•> ' f i w y ar ** -- ^ **#««** 











with the boundary conditions 


y’ "(a) - Jt 4 y' "(8) - Jt, 

y' (B) - tj 
y (6) ■ * 3 • 

If a < 0, 0 < 0, then (x) > 0 and U(x,p) tends to a solution of 

xy* - y - 0 

with ths boundary condition y( 0) - ij, because the Main Theorea requires the canceling of 
two boundary conditions at the right endpoint and of one boundary condition at the left 
endpoint* But if a > 0, 0 > 0, i.e. b Q (x) < 0, then U(x,p) does not converge at all* 



( 101 ) 

(107) 

( 102 ) 

(103) 

004) 

(108) 
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Ha make the following assumptions] 


1°. x is a real variable. 

2°. p is a real positive parameter. 

3°. The real functions a (x) and b (x) have at least 

v u 

n bounded derivatives in the interval 

a < x < 8 « 


4°. n > m > 0 . 

5°. Bq(x) ■ 1 . 

6°. bg(x) y 0, for all x in a < x < 8 * 

7°. n > X 1 > Ij >•••> \ r > 0 

n > T r+1 > T r+2 >, ‘* > T n > 0 * 

Then the behavior of U(x,p), as p tends to infinity, can be 
found by the following procedure: 

First Step . Find the remainder s of the division of n - m by 4. 
Second Step. Find, in the table on the next page, the values of the 
numbers q and p for the differential equation under consideration. 


© If a - 1, b Q > O' 
or 

a ■ 3# bQ < 0- 

(il) If 3 » 1, b 0 < o' 
or 

s - 3, b Q > 0 

(ill) If s » 0, b 0 > 0^ 
or 

s = 2, b- < 0 y 


n-mt1 n-m-1 

then q » —-—, p - —-— 


.. n-m-1 n-m+1 

then q » —-—, p - —-— 


.« n-m n-ra 

then q » —, p - — 


(113) 


(IV) If s - 0, b Q < 01 


8 - 2, b Q > 0 


n-m-2 n-m-2 

then q » —-—, p - —-— 











■fl 


t 

Mi 

* 



Third Step . 

A) If the differential equation under consideration belongs to 
one of the cases 0 - (iTl) * of the table above, try to cancel p 
of the boundary conditions at the point x “ B and q of the 
boundary conditions at the point x • a, going in each group of 
boundary conditions from those containing higher derivatives to those 
with lower derivatives. This is only possible, if there are enough 
boundary conditions on either side to be canceled. 

B) If the differential equation under consideration is of the 
type IV, proceed first as under A). Prom the remaining boundary 
conditions cancel then those two which contain the highest order of 
differentiation without regard to the endpoint at which they are 
given. It can happen that the boundary conditions to be canceled 

in application of this last rule are not uniquely determined, because 
one would have to decide between two boundary conditions of the same 
order of differentiation. He shall call this last occurrence the "Case 
of Indetermination". 

Convergent Case . If it is possible to apply the rule of the Third 
Step in a uniquely determined fashion, then 

u(x) - lim 0(x, p) 

P*** 

exists and is, in a < x < 8, that solution of the differential 


equation 


M(y) -0 


which satisfies the boundary conditions not canceled in the Third Step 
of this rule, provided the following two conditions are satisfied! 


He use the circles around these numbers 
distinguish this division into four cases 
introduced presently. 


, writing 0, (il) , Qy) , (w) , in ord 
from another division into two cases X, 


order to 
X, XX to be 


**~*er?~ 







8°. If the boundary conditions not canceled in tha Third Btap 
of this rula ara raplacad by tha corresponding homogeneous 
boundary conditions, than tha problem determined by thaaa 
boundary conditions and tha diffarantial aquation M<y> • 0 
haa only tha solution y( x) =0. 

8°. Ho two of tha boundary conditiona cancalad in tha Third 
Stap at x « a have ordara of diffarantiation that ara 
congruent modulo n - a, and tha saaia is true for tha 
boundary conditions at tha right and point. 

Divergent Casas . If tha rula of tha Wiird Stap cannot be applied 
because at one of the endpoints there ara not enough boundary conditions 
to be canceled, than U(x,p) will, in general, not converge, as p ♦ «. 
The proof for tha divergence given in this investigation la valid only 
under two assumptions analogous to 8° and 9°, which for their formula 
tion require an additional remark: 

Fourth Step . If the rule of the Third Step cannot be applied because 
the number of the boundary conditions at one endpoint is smaller than 
the boundary conditions that would have to be canceled, then cancel all 
the boundary conditiona on this side and so many boundary conditions on 
the other side (going, as before, from higher to lower order of 
differentiation) that m uncanceled boundary conditiona remain. Then 
we make the assumptions: 

8°'. If the boundary conditions not cancelled in the fourth step 
of this rule are replaced by the corresponding homogeneous boundary con¬ 
ditions, then the problem determined by these boundary conditions and the 
differential equation M(y) - 0 has only the solution y(x> = 0. 

9 . Ho two of the boundary conditions canceled in the fourth step 

at x “ a have orders of dlffsrantlation that are congruent modulo n - m, 
and the same is true for the boundary conditions at the right endpoint. 
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T 


'MKfil 


i 

£ 


Conclusion In the divergent case. 


1) If the rule of the third step cannot be applied because at one of 

the endpoints there are not enough boundary conditions to be canceled, and 

If conditions 8° and 9° are satisfied, then 

lim U(x,p) - ±«, for all x in a < x < 6 . 

P-*- 

2) If the rule of the third step cannot be applied because of indeter¬ 
mination, and If 8° and 9° are satisfied for each of the two possible 

ways of applying the cancellation rule, then U(x,p) does not converge, as 
p ♦ - but rsixains bounded. 


4, The reader is advised to check the examples given in fl in the light of the Main 
Theorem, in example 2, in particular, we discussed a case in which assumption 8° was not 
satisfied. We now give an example where assumption 9° is not satisfied: 


Example 4. n ■ 5, m « 2, b n > 0. 


y* " (a) - t. 


y* (a) - l A 


y I4> (8> - t. 


y‘ <&> - t. 


y <o> - * 5 

Here n-m ■ 3, hence s - 3. From table (113) we find q ■ 1, p « 2. 

The two boundary conditions that are to be canceled at x » 8 have the orders of 
differentiation 4 and 1. But 4=1 (mod n-m), in this case. This means, assumption 9° is 
not satisfied, and the Main Theorem does not apply. However, if b Q < 0, then 9° is 
satisfied, and we can be sure of the convergence of U(x,p). 

It is an open question whether U(x,p) can converge even if 9° is not satisfied. It 
seems unlikely to the author that the Main Theorem remains valid in those cases. 

The next example is of the type which we have called the case of indetermination. 



f *• **!f-*F- 











Example 5 . n - 3, ra - 1, b Q > 0, 

y* ( a) - i 2 y<e) - t, 

y («) * * 3 

Here, n-m » 2 and therefore 8-2. Table (113) shows that this is the case , and 
that p - q - 0. The rules of the Main Theorem would require the canceling of the two 
boundary conditions Involving the highest derivatives. This cannot be done in a uniquely 
determined way, since y(0) - just as well as y{ a) - might be canceled In addition 

to y'(a) - . Hence, this is the case of indetermination, and U(x,p) does not 

converge. 

The rest of this chapter is devoted to the proof of the Main Theorem. 

§3. Asymptotic Solution of the Differential Equation L(x,p) - 0 

5. As pointed out in the introduction the principal tool of our proof of the Main Theorem 
is the theory of asymptotic solution of differential equations involving a parameter. He 
begin by defining what we shall understand by asymptotic equality in this investigation. 
Definition : The functions f(x,p) and g(x,p> are said to be asymptotically equal in an 
interval a < x < 0, if 


.. E(x,p) 

f (x, p) - g(x, p) + -. 

P 

Here a > 0, (but not necessarily an integer), and E(x,p) is a function such that there 
is a positive real number R so that |E(x,p)| is uniformly bounded for o < x < 0, and 
p > R. 

If a function f(x,p) is asymptotically equal to a function F(x) independent of 
p, we shall write 

f(x,p) - (F(x)J 
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1 


Note that the symbol (F(x)J does not describe the function f(x,P) uniquely, 
correct to conclude from 

f^x.P) - [F(x)J 

and 

f 2 (*,P> - lr(x» 


that 


fj(x,P) - f 2 (x,p) . 


It is not 


6. Using Noaillon's method the following theorem can be proved. 

Theorem 1 ; If the assumptions 1° - 6° of the Main Theorem are satisfied, then the 
differential equation (101) admits a complete set of n linearly independent solutions of 
the form 


O u (*,P) 


f 




tn(x)]. 




(V 


1,2,...,n-m) 


(u(x) ] , (V- n-m+1, n-m+2.....n) 

^ v-n+m 


(114) 


(115) 


Here we are using the following abbreviations 

i/ n __ 

1 ) a — | p | 

2) ♦,(x), ♦ (x),...,4 <x) 

I i n-ra 

n-m values of 


<-b Q (x)) 


1 / 


n-m 


arranged in such a way that 


Re(^) > Re<$ 2 > • 


3) 


a (Ob (E) - b (t) 

■r a ' 


■) 


n(x) 


(116) 


(117) 


(118) 
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■*■-■**!•**? 


4 

4 i 

* I 

fc*l 

* 


4) The function* u^(x), (u - 1,2,...,■) are any n aolutiona of the differential 


equation 


M(y) - 0 


forming a complete linearly independent system of such solutions. 

The equations (114) and (115) may be formally differentiated at least n-1 times. 


n^ i> Cx,p) 


for i - 0,1,...,n-1. 




l* v (x)*n(x)] , (v - 1,2,. •. ,n-m) 


* u v-n+m* ' * v " n-m+1 » n-sM-2,...,n) 


7. Remarks. 


As roots of one and the same real function, the complex functions *Mx) are of a 
particularly simple structure. If b Q (x) < 0, then the *Mx> ere obtained by 
multiplying the (n-m)-th roots of unity by the factor 


A similar relation holds when b 0 (x) > 0. More precisely: Set 


, if b Q < 0 


1 / n _ »1 

lb Q (x) | e n-m , if b Q > 0 


Then the «^(x) are, - in different order, - equal to 

. , .2 , n-m-1 

k, ke, ke , ...» , kt 












Furthermore, sine* the * y (x) art the root* of e reel nuaber, the sequence of the ^(x) 
conalste of pelre of conjugate nuabera. Hence, In (117) the •>" end the signs 

alternate. 


§4. Outline of the Asymptotic Solution of the 
Boundary Value Problem 

8. He try to repreeent the eolutlon of the boundary value problem defined by (101) and 
(108) in the fora 


n 

O(x.p) - l c (p) 0 (x,p) . (125) 

vl v v 

If auch a eolutlon exists, then the function c y (p) are solutions of the system of linear 
equations 

n 

l c v (p).L 1 (O y ) - t t , (1 - 1,2,...,n) . (126) 

Using (108) and substituting for the U(x,p) the expressions (114) we find, for 
v * 1,2,...,n-a. 


w - 


X i °"u X i 

a e (8)n(8)J. (i - 1,2,...,r) 

^ ^ , (v ■ 1,2,...,n-m) 

o 1 lF v X (<*)n(«)l. (1 * r*1, r+2,...,n) 


(127) 


where 


»„ - /f *U>dC . (v - 1,2,...,n-m) . 

V d V 


For v « n-mM, 


,n we have, because of (115), 


W ■ t Wn- ,i 


(128) 


( 128 ) 











From the inequalities (117) it follows that similar inequalities hold for the 
quantities w^ (v“ 1,2,... ,n-m), i.e. 

Re(w„) > Re(w ) >•••> Re(w ) . (130) 

i i n-m 

In order to find the c (p) of (125) from (126) we have to calculate the determinant 


L^U,) 

L^Oj) 

• • • ♦ 

LjfO,) 

• • • 

l 2 (u 2 ) 

• • • 

• • • • L 2 (U n > 

• • • 

w 

• • • 

w 

.... L n (U„) 


and the determinants A^tp), (v “ 1,2, 


,n), obtained by replacing the v-th column of 


A(p) by the column 


The coefficients c y (p) are then given by 


A (p) 

c (p) “ . y ■ ■ 

v A( p) 


{v « 1,2,...,n) 


(1.32) 


and, if we substitute (132) in (125) we obtain the form 


n A (p> 

0(x ' p) ” \ MW V X,p) 


for the solution U(x,p) of our boundary value problem. 
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9. Our aim is now to calculate, with the help of tha asymptotic expressions (114), (IIS) 
and (127), (129) the asymptotic value of the right number of (133). The first and most 
important part of that calculation consists in finding the asymptotic value of the 
determinant dtp). The asymptotic calculation of the d^( p) does not offer new 
difficulties. 


$5. The Asymptotic Value of Mo) 

10. As a consequence of (127) and (129) all the terms of the expansion of dtp) are 
obviously of the form 


1*1 d 


S dV 


with real S and real or complex V and K. 
Definition: Two expressions 


K, - t»t 1 l a 


S, dV, 


S dv 

K 2 - (k 2 l a e 


with 


I* 0, k 2 ? 0 

and real S, and s 2 will be said to be of equal order of magnitude if 

Retv,) - Re(Vj) 

and 


s, - S 2 


If 


Re(V,) > Re(Vj) 


Re(V,) - Re(V 2 ), but 3, > Sj 

then K, is said to be of higher order of magnitude than K 2 , and vice versa. 

If R, is of higher order of magnitude than Kj, then we can obviously write 

s dv 

K 1 ♦ k 2 ■ ® • 
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I - 


The sum 

- V ,, 

l [kj o S • 0 34) 

(I) 

of all the t mtmM of highest c ler in A(p) is the asymptotic expression of A(p) for 

large p, unless all the V are alike and 

Hi 

Ik - 0 . 

L m 

(i) 

In this latter case (134) reduced to 

. . s v i 
(0] o • 

and an asymptotic calculation of &{ p) would have to take into account terras of lower 
order of magnitude in the expansion of Mp), as well as the later terms in the asymptotic 
solutions of our differential equation. The exclusion of this exceptional case from our 
theory will compel us to Introduce the conditions 8* and 9* of the Main Theorem. 


11. If the values (127) and (129), for the (U are substituted into the expression 
(131) and A(p), it is seen that the last m columns of the determinant form the matrix 
of n rows and ra columns 


(L 1 (u 1 )] [L^ (Uj)] • • ♦ [I^lu^)] 

(L^U,)] tt, 2 (Uj)] • • • [Ljtu^)] 


| [Ljjtu^) H.„(u 2 )) • • • [LjjtUjj)) 


( 135 ) 


Ml the minors of this matrix have an order of magnitude not greater than that of 1. 

The elements of the first n-m columns of A( p) are given by (127). In order to 
find the asymptotic value of A( p) we expand A( p) in terms of its n-m first columns 
and investigate the order of magnitude of the minors in this expansion. 
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12. Lemma 1: Let D(p) be that minor of the determinant Mp) which is formed by the 


first n-m columns of A(p, and by its 


th th 

1 1 ' S 


. i 


th 


rows, with 


i, < i 2 < • • • < i n . m 


If then the first h of the numbers i are less than or equal to r, then the order of 


magnitude of D(p) is not greater than that of 


h 

0 I W N 

T \m1 
a e 


(136) 


where 


h n-m 

* - I \ * 1 \ • f 137) 

v* 1 u t#*h+1 p 

Proof i From (127) and (131) we see that each of the first h rows of D(p) contains the 
factor n(B)i while each of the remaining rows contains the factor n(a). Hence, we can 
factor in D(p) the expression 


n-m-h 

n ( S) n (a) 


Furthermore, we see, that the first row of D(p) contains the factor o , the second 

X i X i 

2 

row the factor a , etc., and finally the h-th row the factor o . Similarly, we can 
factor in the remaining rows of D(p) the expressions 

T i T i T i 

h+1 h+2 n-m 

o ,0 , • • • , o , 

respectively. Altogether we can factor in D(p) the expression 


n-m-h 

h, .. , . T 

n (6) n (a) a , 
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where T has the value of (137). Then D(p> can be written in the form 


1 

MS 


X. 

e (8)1 • • 


ow 1 

• • e n (6)1 

n-m 


n-m-h 

D( p) " n ( 8) n («> o 



[*, n ” m (a) 1 


"” m (a) 1 
n-m 


Now we expand the remaining determinant with respect to the minors of its first h rows. 
Each term of this expansion contains an exponential factor, and none of these exponential 
factors is of greater order of magnitude than the one originating from the minor formed by 
the first n-m rows and columns of the determinant of (138). For the exponential factor 


in this term is 


h 

o l * 


and because of (130) no sum of h w 's has a greater real part than £ w^. Hence, the 

v v»1 

order of magnitude of D(p) is not greater than that of 

h 

0 I W v 
T \S>1 
0 e . 

It may be less, for we have to take into account the possibility that 


1 •*“ ' * <fr+ 











V 


h 

o l » 

T u* 1 
D(p) » fO 1 a e 

Remark: if Re(w h+1 ) < Re(w h ), than there la only one term of maximal order of magnitude 
In the expansion of the determinant of (138) with respect to Its first h rows. But If we 
have Re(w h+1 ) - Re(w h ) (compare section 7 and the definition of the w^, formula (128)), 
then 

h h-1 

*«( I - v ) - *•( I «v + Vi) ' 

\^1 ^1 

and we have therefore a second term of maximal order in the expansion of the determinant of 
(140) with respect to its first h rows. But, clearly, these two terms cannot cancel, 
since 

*h * Vi • 

13. The question of finding among all the minors of the first n-m columns of A(p) the 
one of highest order reduces now to the two questions: 

(a) which selection of n-m rows in A(p) leads to an expression (136) of highest 
order, and, 

(b) when does the minor corresponding to this selection actually have the order 
indicated by (136). 

In answering the first question several cases are to be distinguished. For this 

distinction the number of v^ix )'s which have a positive real part plays an essential 

role. Let us call this number p. From the definition of the V> *s it is clear that p 

v 

depends on the sign of bg(x) and on the remainder of the division of n-m by 4. A 
simple calculation, which we omit here shows that p has the values indicated in the table 
(113). 

We distinguish the following cases: 

I. Ret^^) ft 0, for all v " 1,2,...,n-m 

A) r-p < 0 

B) 0 < r-p < m 

C) r-p > ra 
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II. Ral^) - 0, for some v . 

A) r-p < 0 

B) 0 < r-p < m+2 

C) r-p > o+2 . 

Remarks : Remembering the definition of the we see Immediately that case II of (139) 

occurs only when 


n-o = 0 (mod 4) and < 0 


(140) 


n-m = 2 (mod 4) and b. > 0 
0 


Case II Is thus seen to be equivalent with the case (iv) of table (113). In this case 


© 


the 


>1ways exactly two functions * with vanishing real parts. In the 


arrangement of (117) these are 


^(x) and i(x) . 
p+1 p+2 


Casa I of table (139) corresponds to the cases ^i) - of table (113). in these 

cases the number q of table (113) is equal to rt-m-p. It follows, therefore, that the 
condition IC of (139) can be re-wrltten In the form 

(n-r) - q < 0 . 

In case of (113), which we have seen to be equivalent to case II of (139), we see from 

(113) that 


q » n-m-p-2 , 

and therefore lie is also equivalent to (141). Hence, we can say, that in case I as well 
as in case II of (139) 

A) is the case in which the canceling rule of the Main Theorem cannot be 
followed because there are not enough boundary conditions at x ■ 6. 

B) is the case in which the canceling rule can be followed. 

C) is the case in which the canceling rule cannot be followed, because there 
are not enough boundary conditions at x ■ a. 
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14. We treat the case I of (139) first. According as to whether we are in the case IA, 


IB, or 

IC, 

we consider then the 

minor D (p), D. (p), 

a d 

D 0 (p) formed 

by the 

n-m first 

columns 

of 

the determinant A(p) 

and by the rows 




1,2, ee 






(case A),(142a) 

1,2, .. 

.p> 

e e e < 

r+1. 

r+2, . 

e , r+n-*-p e 


(case B),(142b) 


r+2, . 


so, n • 

(case C),(142c) 


respectively. Using lemma 1 we shall show that in each ease the einor thus defined is - in 
general - of the highest possible order of magnitude among all the minors of the first n-m 
columns of A(p). 

Case A) The expression (136) has the highest possible order t if the rows of the minor 

h 

are chosen in such a way that J w has the greatest possible real part. In case IA 

v*1 v 

this means that we have to choose h • r. For h can, - by definition, - not be greater 

than r and, on the other hand, all w^ with v < r have positive real parts in 

consequence of the condition r-P<0* In order to make T in (136) as great as possible we 

n-m 

have then to choose for the remaining n-m-r rows of the minor those for which J T 

ir»h+i u 

is greatest. Since the are arranged in order of decreasing size, this is the case, if 

we choose the rows r+1, r+2, .... n-m. This is exactly what we have done in (142a). 

h 

Case B) Here, taking h - r would not make the real part of J w a maximum 

v-1 

because we have r > p and the sume would therefore include w v '* with negative real 
part. Instead, we have to take h ■ p Including thus all the w^'s with positive real 
part and only those. (137) shows then that taking the first p rows of A(p) gives the 
greatest contribution to the exponent T. For the remaining n-m-p rows we take the rows 
r+1, r+2, ..., r+n-m-p in order to make the second sum of the right number of (137) a 
maximum. This is possible because, in consequence of (139), r+n-m-p < n, in this case. 


"H 1 










Case C) Here we reason as In case B. But since r+n-a-p > n in this case, takinq 
the rows r+1, r+2, ..., n in addition to the rows 1,2, ..., p would not be enough to 
have n-a rows altogether. We oust therefore choose h » r-a, in order to have n-a 
rows, and, as before, we see that T is greatest, if we take the first r-a rows of 
A(p). 


Coaparing (142b) with the third step of the rule of the Main Theorea we see that the 
numbers of (142b) are just the subscripts of the boundary conditions canceled in 
application of the Main Theorem in the cases (7) - ^IIlJ, when the canceling is possible. 
Similarly, (142a) and (142c) contains the subscripts of the boundry conditions canceled in 
the fourth step of the Main Theorem in the cases - ^1^. The reason for this fact 
will appear in the course of our proof. 


IS. The considerations of section 14 are not sufficient to prove that the ainors D^tp), 

0. (p), D (p), respectively, really do have higher order than all the others (compare 

D C 

section 12). To investigate this question let us calculate D (p), D. (p), D (p) 

ADC 

explicitly. Since we are most interested in case IB, which we shall see to be the 

convergent case, we discuss O.(p) first. 

b 

In order to calculate D. (p) we write D. (p) in the fora (138). This means that we 

D D 

Interpret the numbers 

1 1' i 2.^n-m 

of (138) as being the numbers (142b), in the same order. Hence, the number h of (138) is 
in this case equal to p. From the definition of p and the w^ it follows that 

S'fVl 1 < Re{w p> * 

The term originating from the minor of the first p columns in the expansion of the 
determinant (138) in terms of its first p rows is therefore the only term of maximal 
order of magnitude in this expansion. We obtain therefore 

P 

a I w 

n-m-p T. ** v 

D fe (p) ■ (^(8)8 (a) • • B fa J o e (143) 
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with 


and 


! 


x ♦ 


V 


r+n-m-p 

l 


■t 

u 


(144) 


* t (B) • 


• • * ( B) 

P 


*,**( B) • 


• • * MB) 
P 


(145) 




r*"— 

n-a 


v" 1 ** 


„ r+n-a-p, 

Vi ta) 


(146) 


(134) is valid also for p “ n-m and for p “ 0, If M define 
Bjj - 1, for p “ n-a. p(x) is an exponential function, hence 
D. (p) has therefore the order of magnitude of 

D 



Aj, - 1. for p - 0, 
n(B) »• 0, n(u) * 0. 


and 


if A b ^ 0 and 0. 

Lemma 2 : ^ 0 if and only if 

\ i i A^, (mod n-a)> (i,J » 0,1,...,p) (147) 

and Bjj f 0 if and only if 

t i t., (mod n-m)> (lc,i • r+1, ...,r+n-m-p) . (148) 

n 1 


Proof : The statement is trivial, as far as is concerned, if p • 0. For p > 0 we 

remember from $3, section 7, that, with the notations used there, the * v fx) are, in 
different order, equal to 
















I 


k(x), k(x) e, k(x) t , • ... , k(x) e 1 *"""' . (124) 

In tha ordar of (124) the •f^ ara rapraaantad in tha complex plana by a sequence of 
successive points on tha circla of radius |k|. It is aasy to sea that tha numbers 1 ^( 0 ), 
¥> 2 (8)» .•• i v> p < 6) ara than, in different order, equal to tha numbers 

k(0)e t+1 , k<0)e t4,2 , • • • , ktB)®** - * (149) 

where t is a certain integer which is only determined modulo n-m. Substituting the 
expressions (149) into (147) we find 




\^i 

i k(0) 


(tMU. 


(ttl)X 


(t+p) X, 
e 


(ttp) X 
c 


(ISO) 


Now we set 


C v' C ' 

This allows us to write (150) in the form 


(v* 1,2,« ♦» ,p) 


\p1 


^ - * k(6) 


t+1 


t+p 

• • • C- 


t+1 tfp 

5 i ' * • s 


(151) 




where is the Vandermonde determinant of ^Sinca tha 

Vandermonde determinant vanishes if and only if two of its rows are equal, ia zero if 

and only if some of the i; are alike. But from (151) we see that 


means 


C 


1 


or, because of the definition of e in formula (123), 

Xi = Xj , (mod. n-m) 

This proves the part of lemma 2b that is concerned with The proof for is exactly 

analogous, and is therefore left to the reader. 

The reader will readily see that similar results hold for D (p) and D (p). The 

4 c 

only real difference in the reasoning comes from the fact that in these cases there may be 
two terms of maximal order in the expansion of (140). But since these terms cannot cancel, 
this does not essentially affect our argument, we restrict ourselves to stating the 
results in these casesi 


T 0 l W v 

D # (p) “ fft a ) a a . ^ 

r-m 

o £ w 
T ^ v 

D(p) - [Q J a c e 
c c 

where 


r 

I 

X + 

n-m 

l T u 

v*1 

V 

l^r+1 u 

r-m 

l 

X + 

n 

I T 

V*1 

V 

irr+1 “ 


( 153) 


(154) 


(155) 


(156) 


and Q a and Q c are two constants, with respect to which the following two lemmas hold: 
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Lemma 2a i O a ^ 0 If 


and 


£ X^, mod n-m, (i,j « 1,2,...,r) 


(157) 


£ Tj, nod n-m, (k,I - r+1. 


,n-m) 


Lemma 2c: Q_ *< 0 if 
— c 


\ £ \y 

mod n-m. 

<i,j - 1,2,.. 

.., r-m) 

T k * V 

mod n-m. 

(k,i - r+1,., 

>• «n) 


(158) 


Remark : From these lemmas and the rmurk at the end of section 14 we recognize that the 
conditions (147) and (148) are equivalent with the assumptions 9* of the Main Theorem while 
(157) and (1S8) are equivalent with 9*'. Note, however, that our reasoning so far does not 
cover the case II of (139), which we have seen to be the same as the case fiv) of (113). 


16. The cofactor of the minor D^tp) in 4(p) is that minor of the matrix (135) which is 
formed by the m rows of A( p) not contained in D^(p). let 


*b 


Vi'V. 

• * * • Vi'V 

W. 

• • • • L (u ) 
r m 

L r+n-m-p+1* u 1* 

• • • • X* 

r+n-m-p+1 

VV. 

• • • • L (u ) 
n w 

Dtp) is of the form 

D 



(u > 


(159) 


± iy . 

We now introduce the assumption 

^ + 0, in case IB . (160) 

This condition will be seen later to be equivalent with the assumption 8* of the Main 


Theorem in the cases 


Then the term t D b (p> (y is the term of greatest 






















order of magnitude in the expansion of &( p) in terms of its first n-a columns. Hence. 
(143) gives 

T. 0 \ *v 

A(p) “ [Q b • ^1 a . / in case IB / (161) 


where 


n-m-p 

2 b - ±rf(e> • n (a) • • Bjj ? 0 


(162) 


Similarly, if we define and 6 c as the limits - as p of the matrix (135) 

formed with the rows of D(p) not occurring in D (p) and 0 (p), respectively, and if 

B C 

we introduce the assumptions 


6 y 0, in case IA 
a 

6 y 0, in case IB 
c 


then we have, in analogy with (161) 


T 0 l W v 

A( p) - • 6 # ] • o e , 

r-m 

_ o y w 

T V 

A( p) • [Q •6]»o C « V " * 

c c 


in case IA , 


in case IC . 


(163) 

(164) 


(165) 


(166) 


This finishes the asymptotic calculation of Alp) in the case I of (139). 


§6. The Asymptotic Value of the Solution of the Boundary 
Value Problem in the Case IB. 

17. The method used in $5 for the calculation of Alp) can also be applied to the 
determinants A^( p). (For the definition of A y (p), see §4, section 8). 

For v > n-m, A^( p) is distinguished from Alp) only in one of the last m columns. 
In this case all the considerations of |5 remain valid for A y (p)■ if the determinants 
6 , 6 are replaced by the determinants 5^, 4^, 4,^ obtained by replacing the 
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i-(n-m)th colvatn of 4 # , 5^, by the cor responding 


!•••> by 


*n-a*1 *pM 1 r-«+1 

• • • 

* • • 

• ■ • 

• l or • 

r 

• jl • 

rtn-n-p+l 



and analogous formulas hold in the two other cases. 
Hence 


(v > n-m)i case IB) , (167) 


and similarly in the two other cases. 


y p > ry 

& ( P ) ■ ly 


18. For v < nr-m the column of A( p) that must be replaced by the i^'s in order to 
obtain A Ji p) changes the structure of the determinant somewhat. But if we place this 
vth column behind all the others (and change the sign of the determinant, if necessary), 
then we obtain a determinant very similar to A(p)• The only essential difference is 
that a has to be replaced by a' ” ih- 1 (and therefore n-a by n-m-1) and that w^ 
must be omitted from the sequence 

*1<*2' • • * ' w n-a* 
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T 


Again we treat the case IB (the convergent case, as we shall see) first, and since we 

V p) 

are more interested in - — than in A Ji p) itself, we state the result in the following 
font i 

Lemma 3: In case IB, l.e., if 0 < r-p < m, we have 


-X -ow 

o P e U t > v < p 


r+n-m-p , . 

a [ it ]» p < v < n-m , 


where the constants depend on the values of the J^'s and are not necessarily 

different froa zero. 


SI 


1.) v < p. In this case the w^ to be omitted has a positive real part. Hence, the 
formulas for A( p) can be applied to A J p) if we replace 

p by p’ - p-1 
m by m' • m+1 

and the sequence 


by the sequence 

w 'l' w 2.^m-1 

identical with the first one except for the omission of the term since we have 

0 < r-p 1 < m 1 , 

the formula to be used in (161). He obtain 

P’ 

T 0 £ W u> 

V p) " {Q bv * W * 0 bV 6 UPl 








I 


i 


p' r+n-a'-p' 

T m l X * l T 

<i^1 ipr+1 


Note that the reasoning used for the proof of ^ 0 does no longer hold for 

Is a determinant depending on the . It Is defined as the cofactor of the leading minor 

In the expansion of fi (p) with respect to its first n-m-1 columns- 4 is a 

V D V 

determinant of m+1 rows- 

From the definition of p' , a', w' it follows thsn that 


r-.! 


) w* • > w - w 

u • w v 
upI up 1 


p-1 rfn-»*p 

t - y x ♦ j t 

crl irr+1 w 


Comparing these formulas with the expression for M p) we see that 


VI f °bv ‘ *bv ) _X , 

4 “>>" l« b • vJ° 


p . v 


2 .) p < v < n-a. Hare we reason as under 1-), the only difference being that 


p' • p and therefore 


r -• -! 

u* 1 ur i 


p r+n-m-p-1 

T. - ; X ♦ l i 

bu L . u u 

df> 1 ^ r+1 


Hence, In this case 


•> 4 j 

























3 


It is easy to see that the first two sums of the right member of (169) tend to sero, in 
a < x < 0, as p, and therefore a, tend to infinity. Tor, 


R «( / - w ) <0, in a < x < 0, for v < 


f a < 0, in a < x < B> for v > p , 


in consequence of (117), (128) and the definition of p. 

Therefore 

u(x) « lim U(x,p) - £ ^^ u ' (170) 

p+” 1^1 W 

But if we remember the definition of « b , ^ and «Mx), as given in (159), in 

section 17, and in theorem 1 (section 6), respectively, then we see that the right member 
of (170) is just the solution of M(y) - 0 satisfying the boundary conditions not canceled 
in application of the Hain Theorem, and condition (160) is seen to be equivalent with 
assumption 8° of the Main Theorem in the convergent cases (T) - ^II) (table (113)), 


20. Remark : Formula (169) is, in fact, a complete asymptotic solution of the boundary 
value problem. It might be used for a more detailed description of the boundary layer 
phenomenon. It can be written in the following simpler and more symmetric form: 

_T P of mac 

U(x,p) - u(x) + a T £ e (it n(x)J 

\^i v 


i n-ra of * (C)d£ 
♦ o‘ X l e v ( 

\pp+1 


(* v n(x)] 








I 


& • 


“1 

is 


where T « t and X » X are the lowest orders of differentiation occurring in the 

r+n-m-p p 

canceled boundary conditions on each side, An easy consequence of (169) is, e.g., the 
following interesting 

Corollary 1 Under the assueptions of theorem 2 the derivatives 

(t -1) 

U(a,P), U* (<*,P), < • .,U(a,p) r+n_ “"P 

converges to the value of the corresponding derivatives of u(x) at x “ a. The next 
derivative, 

n A __ ) 

r+n-m-p 


0 (a,p) 


(t > 

r+n-m-p 


converges, as P + •, but in general the limit is not equal to u (a) . All the 

higher derivatives of U(x,P) tend to infinity at x ■ a. An analogous statement can be 

(X J 

made at x - 0, with u(B,P) as the last convergent derivative. 

In less precise language we may express the statement of this corollary by saying that 
the last canceled boundary condition at each end point determines the derivative of 
U(x,P) in which the boundary layer occurs at that endpoint. 

§7. The Proof of the Divergence in the Cases IA and IC . 

2t. We are now going to show that in the two remaining cases U(x,P) tends to infinity as 

P ♦ To that end it is sufficient to prove that one of the terms in the right member of 

(133) tends to infinity. For each of the first n-m terms of (133) is of the form 

Y a W v (x> 

K v a e 


Y 0 Wjtx) 

If one of these terms, say K^o e tends to Infinity, the whole sum can remain 

hounded only if the sum of all terms of the same order of magnitude as this particular term 
vanishes identically. This would require that at least one other term has the same 
exponential factor as Kj. We would have, e.g.. 
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I 


a 


! 


w^(x) » Wj(x) . 


Now all tha H (x) are of the fora 


V*> ■ V v* /a V £)dC 


where the constant V originates from the factor . while ( x <f (E)d£ is the 

v A (p) J a v 

contribution of U y (x,p). But an equation like 


V 1 + /a V C)d5 “ v 2 + * 2 <5)d5 


is impossible, since no two <f (x) are equal. 


22. Proof of the divergence in the case IA: 

Let us calculate the determinant A (p). In this case the formulas for A(p) can 

r +1 

be applied, if we replace 

p by p' ■ p-1 (since r+1 < p) , 
m by m' - m+1 , 

and the sequence 


w 1' w 2. w n-m 


by the sequence 


w 1' w 2'. w n-m-t 


obtained by omitting the term w f+ ^ from the sequence of the w^. As we have 


r-p' < m 1 


the formula to be used is (165). we obtain 


j u 

. . . . , a,r+1 ufI 

Vi 10 ’ ■ tQ a,r+1 * 4 a,r+1 5 ° 








It follows that 


r n-m' 

T . , “ [ X + l T 

4ir+1 «o w 11 


r 

I 

Uf*1 


r 

l 

w*i 


and 


a,r+1 


r n-m-1 

l X * It 

*1 “ r*n " 


For 


W» 

A( p) 


we obtain 


r+1 


A ( 


<p) fta.r+l 6 a,r+l1 ”V« 

<» "I fi . 5 . J° 


and the r+1 at term of (133) becomes therefore 


w> 

A (p) 




“a, r+1 


a.r+1 


nC it) 


- T 


r+1 


(0«C 


Since r+1 < p, the exponential factor of the right member tends to Infinity as p ♦ «*• 

If we can prove that the expression in brackets does not vanish, then the divergence of 
U(x,p) is assured, n(x) does not vanish in a < x < 6 (see (11S))« For fi 4)r+1 we 

can prove that the following lsmaa is true. 

Lemma i t Q # r+1 ? 0, if assumption 9*' of tho Main Theorem is satisfied. 

Proof : An almost literal repetition of the arguments of $5, applied to A^^(p) instead 
of to A(p), which we shall omit, shows that is different from sero, if the 

following determinants do not vanish: 
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X. 

*,( 8 ) 


and 


a,r+1 



* 



( 6 ) 


a,r+1 


Vi Vi 

* A a). * (a) 

r+2 n-m 

T n-*-1 Vm-1 

* ,<a). * (a) 

r+2 n-m 


To these determinants th« reasoning of luma 2b can be applied, leading 

6 ., Is a determinant one column of which is formed of m+1 of 

• # rri 

1,. 8 .. will be cero for certain exceptional vlues of these t.. 

assumption that the have these exceptional values would not be sufficient to guarantee 
the convergence of U(x,p), for the r+1 at term of (133) will, in general, not be the 
only one that has an exponential factor tending to infinity* We shall omit the not 
difficult problem of finding sufficient conditions for the convergence. We will assume 
instead, without mentioning it each time, that the do not have these exceptional 

values. 

From the preceding considerations it follows that U(x,p) is in fact divergent in 
case IA, provided that the assumptions of ths Main Theorem are satisfied. 


exactly to lemma 4. 
the n numbers 
But even the 


23. The proof for the divergence of U(x,p) in the case ZC is analogous to the proof for 
the case ZA, if the (r-m)-th term of the sum in (133) is considered, instead of the 
(r+11st, on which the proof in the case ZA was based. One obtains ths expression 
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A r-, (P> 

TIFT 


-P 


-S^ 


* n(j. )J 


-A _ 0{P * 


(E)d£ - w } 
ot r-m r - m 


c c 

from which the divergence of U(x,P) follows as in the case IX, since Re(^ r _ m ) < 0. 


$8. The Case II . 

24. The case II, i.e. the case when Reti^p^) “ Re ^p+2^ ” °» requires a special 

discussion, because in our reasoning in case I we assumed repeatedly - notably in section 

14 - that every 1*^(0^), with i < r, is either of greater or of lower order of magnitude 

than all L^lUy) with i > r. But this is no longer true in case II, for the L t< D p+ i) 

°w i ow 

and I‘ 1 (Op +2 >» because the exponentials e and e ™ have always the absolute 

value 1. We must therefore modify our considerations for case II, from section 13 onward. 

We remind the reader of the remark made in section 13, to the effect that case II of 
(139) is equivalent with case (iv) of (113). 


25. The case IIX. 

Re-reading sections 14-22 one sees that no modification of the proofs for case IX is 

necessary to obtain the proofs for case IIX. Xll arguments remain literally the same. The 

Ow i «w 

reason is that, for r < p, e ** and e *** do not appear at all in the asymptotic 
expressions for A(p) or (P). 


26. The case IIB. 

This is no longer true for case IIB, defined by the inequalities 0 5 r-p < m+2. In 

order to find the minor of maximal order of magnitude among the minors of the first n-m 

columns of A(p) we go back to lemma 1, section 12. Whether the number h in (136) is 

h 

chosen equal to p, p+1, p+2, does not have any influence on the real part of I w , 

h v-1 

but any other value of h leads to a \ w with a smaller real part. Which of these 

v»1 

three values has to be chosen for h in order to make the order of (136) a maximum depends 

therefore on (137). h must be that one of the three numbers p, p+1, p+2 for which the 
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X^ and the can b« chosen to w to oaks T a maximum. Let ua call this maximum 

v u 

T^ 11 *. As in section 14 it is clear that nust belong to T^ 11 ^, as well 

M T r+1' T r*2.V^p-2‘ 

With respect to the two retraining tents to be chosen we can only say that they nust be 
the two largest of the remaining X*s and tJs. This proves that the minor of greatest 
order of nagnitude among all minors of the first n-a columns of &( p) is in this case 
forned by the rows 

1,2,. . .,pi r+1, r+2,.. ..,r+n-n-p-2 , (172) 

and two nore rows, which must be those corresponding to the boundary conditions containing 
the highest order of differentiation excluding those rows already contained in the sequence 
(172). A comparison with part B of the third step of the Main Theorem shows that these are 
just the rows of h[ p) belonging to boundary conditions that must be canceled in 
application of that rule. We know from the proof in case I why this is sot the rows of 
A( p) appearing in the minor of maximal order ot the first n-m columns of M p) are 
just those that do not occur in the cofactor of this minor, and this cofactor determined 
the boundary conditions satisfied by the limit of U (x, p), if it exists. 

The exceptional case of indetermination occus, when the two additional rows after the 
rows (172) have been chosen, are not uniquely determined. We shall assume, for the 
present, that we have to do with the regular case. The case of indetermination will be 
treated in J9, 

Our rule is so formulated that it takes also care of the case that X .. or X ., do 

pt i p* . 

not exist because r«p or r « p+1. 


27. The minor of maximal order or magnitude among the minors of the first 

D 

n-m columns of &(p) can now be written 


.(II) 


d £ II ’<0> 




H 

° i 

V»1 


(173) 
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where 


,<II) 


H r+n-m-H 

I v i 

w*1 ivrM 


(174) 


and H is equal to p, pH, or p+2, as the case may be. 

As in case I we can then prove that Q* 11 * ? 0, provided assumption 9* is satisfied 

D 

The cofactor of D* 11 ^(p) will again be of the fora 
D 


iC 1 * 1 


where A^ 11 * is defined in analogy with Assuaption 8* of the Mein Theorem assures us 

again that 


* 0 • 


Finally, we find for A(p), similar to (161) 


A(p) - c b 


(id o i » 

- Vi 


in case IIB 


(175) 


28. As in case I, we could now calculate the determinants A^{p) by the method used for 
Mp). But since we are only Interested in proving that 


V p> 

A (p) 


« v (*.p) 


tends to zero, for v < n-m, we shall not calculate the analogue of (169) for this case. 

As in proof of lemma 3 we consider first the A (p) for v < p. For these v the 

°"v 

determinant A (p) does not contain the column with e of A (p). It follows from 

v V 

lemma 1, that then the minor of maximal order among the minors of the first n-m-1 columns 

of A (p) cannot contain as factor an exponential of higher order than 
v 








"~T 


X 


a 7 w - w 

A « v 

• ap1 

This is thsrsfors also ths sxponsntial of maximal order that can occur in the asymptotic 
expression for A (x). Hence/ 


has the exponential factor e (and possibly exponential factors of order 1)/ and 


T7p7 °v l *' pl ' U<p) 


has an exponential factor whose exponent is 


< »<«>««-v 


an expression which tends to sero as 0 ♦ «■ 

V p> 

If v > p*2, we prove eimilarly that does not contain any exponential factors 

(except possibly exponentials of the order or magnitude of 1). Hence the asymptotic 
expression for 


y p) 

^H v (x,p) , ( v > pf 2) 


contains an exponential factor whose exponent is 


o * v U)d5 / ( v > p*-2) , 


and this expression tends to sero as o ♦ 
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fcS 


29. The two remaining terse of (133) require a more careful analysia. We treat only the 
case v - pH, the caae v « pt2 being almost identical. Similar as in section 26 we ask 
which choice of n-m-1 rows of A_, . (p) leads to the minor of highest order among the 

ptl 

minors of the first n-®-1 columns of A_^.(p). This minor will have either 

ptl 

a\ w a 7 w + w 

« <>> .«* P+2 

<tf*1 IV*1 

e or e 

as exponential factor. In both cases the asymptotic expression for 

A < p) 

AM~ Vl tx ' p) 

will contain an exponential factor of the order of 1. But the asymptotic expression for 

T' 

A ( p) has also a factor of the form a , and similarly as in section 26 we conclude 
ptl 

that T' is a sum which contains the terms X 1 , X^,•••» X^ and the terms T f+1 , 

t To these terms one more term has to be added (not two terms as in 

r+2 rtn-m-p-2 

the case of A( p), because A , (p) must be expanded in terms of its first n-m-1, not 

pt 1 

n-m, columns). This term must be the largest of the remaining XJs and tjs. Finally, 

a , i _(H> 

A <p) T '-T 

in the asymptotic calculation of \ . , we have to form g , and this will be 

Vi (p) 

the power of o occurring in the asymptotic expression for . , V J4 (x,p). This 

A (p) A(P) ^ 

proves that the order of magnitude of \ . U . «(x,p) is not greater them that of 

Alp) p+1 


where s is the smaller of the two terms chosen for T. 7 , after X.. *X*,.# •,\ and 
——— d i 2 p 


ww 


, t _ have been selected, s may be any of the numbers 

n- xsr p“ 2 

V 1 ' V * 2 ' T r+n-ra-p-1 ' T r+n-m-p 


In order to prove that the ptl st term of (133) tends to zero it remains only to show 
that none of the numbers (176) can be zero: 











*> ■ “ X . In thi» case the other on* of th* two last numbers chosen for T* 11 ' 

P+1 b 

■uat be T r+n _ m .p_ 1 > because. If It were X^. e would not be the smaller of the two, 

and it cannot be x since 

r+n-*-p 

T r+n-m-p < T r+n-m-p-1 


and the two chosen numbers must be the largest of the numbers (176). Since s must be the 
smaller of these two numbers, we conclude that 


On the other hand. 


X < x . , 

pH r+n-m-p-1 


V+1 > T r+n-m-p 


I< X “ 0. X , must be the last of the numbers X, hence 
P+* p+■ 1 

p+1 “ r . 

It then follows from (177) that x cannot exist, in other words 

r+n-m-p 

x . . « x 

r+n-m-p-1 n 


(178) and (179) imply 


r+n-m-p-1 - n 


m ■ 0 , 


a case excluded from our considerations. Hence, X . . > 0. 

p+1 

b) a « X^,. Similarly as in a) it follows that 
P+2 

^p+2 > T r+n-m-p-1 

X ., ” 0 would imply p+2 - r and r+n-m-p-2 » n, hence 
P+* 

ra ■ 0 . 


c) and d) a - x , or s " x . .We leave the proofs in these cases to 

rTn“m*p'* 1 r+n*B“p 

the reader, since they follow exactly the pattern of the proofs in a) and b) with the same 
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30. Thus we have finished the proof of the following theorem. 
Theorem 3 i If the assumptions 1* - 9* of section 3 are satisfied, if 

0 < r-p < mv 2 , 


and if 


n-m s 0 (nod 4), and b^ < 0 
or 

n-m = 2 (mod 4), and b Q > 0 , 

then, as p ♦ «, the solution U(x,p) of the problem defined by (101) and (108) converges 
in a < x < 8 to that solution of the limiting differential equation M(y) - 0 which 
satisfies the m boundary conditions not canceled in application of the rule of the Main 
Theorem for the case (xv), unless we have to do with the case of indetermination. 


31. The case IIC. In case XIC (i.e., if r-p > mt2) one can again, as in section 23, 
consider the (r-m)th term of the right lumber of (133) and prove that its asymptotic 
expression contains the exponential factor 

<>{/%< Od5 - w ) 

a r-w r-m 

0 • 

This is sufficient to show that this term must tend to infinity, since this exponential 
tends to infinity more strongly than any potter of a may tend to zero. 


§9. The Case of Indetermlnatlon . 

32. let us assume that we cannot decide in a uniquely determined way, which is the last 
row of A( p) to be chosen in order to obtain the minor of greatest order of magnitude 
among the minors of the first n-m columns of &(p). This is the case which we have left 
aside in section 26. We have seen there that this occurrence means that we have to do with 
the case of indetermination of the Main Theorem. 

In this case we must have 

X t“ \ 

where i is one of the numbers pM, p+2 and k one of the numbers r+n-m-p-1, r+n-m-p. 
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Furthermore, X^ and must ba tha aacond and third in aisa of the nunbera and 

■t., excluding tha numbers X.,X_,...,X and t ., _. For simplicity 

i id p r+i r +2 zrn*Bhp "2 

lat ua assume that t * pM. {Tha raaaoning ia the same for t - pf 2.) Denote by D 1 (p) 
and ( 0 ) tha ainora of the firat n-ra columns of tha determinant Atp), with tha rows 

1.2.. . .,p, p+1» r+1, r+2,...,k-l 

1.2.. ..,p; r+1, r4-2,...,k-1, k , 

raapactivaly. Than D^tp) and (p) are both minora of the same maximal order, and we 
cannot raaaon aa in aection 16. But, uaing formula (138) we can write 


ow 


D^p) * F( p) e 


P*1 


[d,] 


D 2 (p) - Tip) 


td 2 ) , 


where 

^ * v + I t 0! w 

Ft p) - O'" 1 ^ “a - 1 

d 1 and d 2 are certain conatanta, which can eaaily be proved to be different from zero by 
the method of aection IS, provided assumption 9* is satisfied for each of the two ways of 
canceling corresponding to D^tp) and to 0 2 (p), respectively. 

The cofactors of the minors D^ {p) and D 2 <p) in At p) can be written in the form 
[5j] and [Aj], where 6^ and A 2 are certain constants. If assumption 8* is 
satisfied for both waya of applying the cancellation rule of the Main Theorem, then we can 
be sure that 

A, f 0, A 2 + 0 . 

expanding At p) in terms of its first n-m columns we find then 


At p) -D^pMA,! tD 2 tp)lA 2 ] - Ft p) {e d^] ± ld 2 A 2 l} . (180) 
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33. While in all cases so far the last a terms of the sum (133) were convergent, 
the convergence or non-convergence of U(x,p) depended entirely upon the first n- 
of (133), in this case the last m terms of (133) will, in general, be divergent. 
In fact, let n-m < j < n. Then we find immediately in analogy with (100), 

oh 

4j( P ) *F(p){e ^’td^j) + (d 2 «2j n ' 

where [6,^] and (5^1 are the determinants obtained instead of (6^) and S 2 ), 
j-th column of A( p) is replaced by 


t 


(180), (181) and (115) show that the j-th (J n-m) term of (133) is equal to 

ow 


Mp> 


d(rt v x ' p) 


>1 


Of 


l Vlj 5 + [ V21> . . . . 

-*- J- [u (x) ] 


P *' 1 [d, 6,] + [d 2 S 2 ) 3 " n " B 


and that 


n A. ( p) 




OW n n 

« l [d t « u (x) ) + l [A 2 & 2 u (x) 1 

i**n-m _ 3 j-n»ra j-n+a _ 3 j-n-m 


0*» . 

e ^ Id, 6,1 


Ed 2 5 2 l 


Because of the oscillating factor e pf1 this expression can converge only if the 
determinant 


i d , a ,j u <*> 

j-n-m J j-n-o 


Vi 


l d 2 « 2 u (x) 
j-n-m 3 j*n+m 


V 2 


vanishes, i.e. if 


so that 
m terms 

(181) 

if the 


( 102 ) 

(183) 


- 49 - 








n 4 2j | 

l u (x) - 0 . 

j-n-m « 1 i 2 [ j-n+m 

Sine* the u y (x), (u - 1.2,...,m) are, by assumption, linearly independent, thie 
leads to the n conditions 

$ 2 5 ij “ 5 ^2j “ ° ^ • (184) 

The left sides of (184) depend on the numbers t,, , (184) represents therefore a 

12 n 

set of conditions on the prescribed boundary valuese As in the previous cases we assune 
that the tjs do not have the very special values required by (184). Then, the value of 
(183) oscillates, as a * •, but remains bounded. 

As to the first n-m terms of the sum in (133), it can be proved exactly as in the 
regular case IIB (section 27-29) that they all tend to zero. 

This completes the proof of the non-converqence of U(x,p) in the case of 
indetermination and also the proof of the whole Main Theorem. 


34. Remarki Going over our whole proof we see that the assumption m > 0 was not used at 


all in the proof for case I, i.e., for the cases 


©- © 


of (113). Hence, all our 


In case II, i.e. (IVJ, t* 1 ® hypothesis m > 0 was used only to exclude a rather 


results in these cases remain valid for m ■ 0 also. 


special occurrence in section 29. it would not be difficult to formulate a general theorem 


for the case m - 0 also. 
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FURTHER RESULTS IN THE CASE n-m - 1 
$1. Generalization of the Boundary Conditions . 

35. In the special case n-m - 1 It is not difficult to replace the boundary conditions 
(108) by the more general ones 




n (v-1) n < w-1 > 

I o iv y (a) ♦ I 8 y (8) 

v-1 u-1 


(1 m 1/2/... ,n) 


( 201 ) 


where the constants ci^, 8^ and are only restricted by the condition that the n 

boundary conditions are independent and compatible. 

We note first: If, in (201), the L,(y. _. .._ . _ _ . . . . . .. 

i J ) and the are subjected to the same 

linear transformation with constant coefficients and non-vaniBhing determinant, then the 
resulting equations 


n (v- 1 ) n (vi - 1) 

LMy) - l a' y (a) + J B* y (0) - «j (202) 

v-1 v-1 W 

constitute a set of boundary conditions equivalent to (201) in the sense that a function 
satisfying (201) satisfies also (202) and vice versa. 

From this remark we see that we can assume without loss of generality that not all the 
or all the are zero, because in that case the boundary conditions (201) would 

be equivalent to ordinary initial conditions, which are a special case of boundary 
conditions of the type (108). 

36. Denote by s the greatest value of v for which at least one of the a, is not 

lv 

zero, and by t the greatest value of u for which at least one of the B^ is not zero. 
There exists a non-degenerate linear transformation with constant coefficients transforming 
(201) into the system of equivalent boundary conditions 
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(203) 


o*i 


n 

■ I 

vl 


a. 

a i v 


< v-1> 
y (a) 


l 

ifl 


a ip 


C M-1» 
y (8) - 


A 


having th« property that 

a. i* 0. cu^ ■••••* a “0 . 

a It r ait a ne 

Such a transformation can be chosen in many ways. Similarly, there is a transformation 
changing (201) into 


6 


n (\>-1) n <p-1) 

*1 S J, 0°iv y (a) + J, 0 6 iy y (6> 


)*■ 1 


0 l i 


(204) 


having the property that 

S 6 1t * °' 0 0 2t B K nt 

For n-m *• 1 aquation (114) and (IIS) reduce to 


0 . 


-p £y° d * 

U^x.p) - e (n(x)l 


U (x.p) - [uU>, , (p- 1,2,...,n-1) . 


Let us assume first that b 0 (x) > 0. Then we use (203) instead of (201) and 


immediately that 


yiy - ± p*" 1 (b Q (o0 1 • n(o>] 


is of greater order of magnitude than all the other L^(U^), (i “ 2,3,...,n). Solving 
equations (125) asymptotically in this case, by the method used in the first chapter, we 
see that we can reason exactly as there. Of the assumptions of the Main Theorem we need 
only 1*, 2*, 3*, 5*, 6* and an assumption corresponding to 8* which states that 








aV u 1>' 


' _ L 2<“ n > 
a 2 n 


L (u ), 
an 1 


* • > 
an n 


+ 0 


(205) 


We see then easily that the solution of (101) and (201) tends, for n-m * 1, and bp > 0, 

to that solution of M(y) ■ 0 which satisfies the boundary conditions 

L,(y) - t. , (1 - 2,3,...,n) (206) 

a 1 a 1 


but not the boundary condition L (y) - 1. except for special values of the 

a l ex i 

a similar reasoning in the case bg(x) < 0 we find that in that case u(x) ■ 
satisfies M(y) - 0 and the boundary conditions 


i’s. Using 
lira U(x, p) 

P-m. 


0 V y) ■ 8*i • (i “ 2 ' 3 


(207) 


provided 


V 


B L 2 ( u 1 )' * * * ' bW 


e L n (U 1 ) ' • • • ' bW 


+ 0 


(208) 


This result can be formulated in a somewhat more symmetrical form. TO that end note 
that the boundary conditions (206) do not involve any more the highest derivative at x ■ a 
occurring in (201). any linear combination of the equations (201) which does not contain 
(s-1) 

y (a) must be linearly dependent on equations (206) and no linear combination of 

(s-1) (t-1) 

equations (206) contains y (a) • Similarly for (207) with respect to y (8) • Hence, we 

can state the following theorem. 
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Theorem 4 » if n-m “ 1 and if tha conditions 1», 2*, 3*, 5*, 6* of tha Ha in 
Thaorea of chaptar I as wall aa (205) and (208) ara aatiaflad, than tha function 
U(x,p) satisfying (101) and tha boundary conditions (201) tanda with incraaaing 
p to a solution u(x) of *K y) “ 0. According as bg > 0 or b^ < 0 
tha function u(x) satisfies all boundary conditions that dapand linearly on (201) 
and do not contain tha highest derivative at x - a or x “ 0, respectively, 
occurring in (201). 

Remark ; The conditions (205) and (208) can be formulated in a way independent of tha 
particular choice of the fundamental systaa u^(x), (u • 1,2,...,n-1), by saying that we 
assume that only the function u(x) =0 satisfies tha differential equation M(y) « 0 and 
the homogeneous boundary conditions corresponding to (206) or (207), respectively. 

Example: n - 3, a - 2, 

y( a) - y' ( a) + y* • ( 8) - i, 

y(n) + y* (a) + y(0) - y* (0) - * 2 

y«(a) - 2y' ’ ( 8) - t 3 

If b Q > 0, then lira U(x, p) satisfies the boundary conditions 
P*“ 

2y( a) + y( 0) - y* ( 8) ♦ y* ’ ( 0) - i, + t 2 

y(a) + y(0) - y’{0) + 2y"(8> - » 2 " » 3 • 

But if b Q < 0, then lim U(x,p) satisfies the boundary conditions 

p+- 

y( a) ♦ y* (a) ♦ y( 8) - y’ ( 8) - * 2 
2y(a> - y'(a) - 2t, + t 3 
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$2. The "Stretching* of the Boundary Um 
37. For the relatively simple types of boundary layer problems with which this 
investigation is concerned we have been able to develop a method that allows us to 
calculate asymptotic expressions for the solution of the boundary value problem (compare, 
e.g., formula (171)). From these asymptotic expressions one can easily obtain all desired 
information about the behavior of the solution of the boundary value problem near the 
endpoints for large values of p. (Compare, e.g., the corollary in section 20). 

In the more complicated boundary layer problems occurring in physics such complete 
asymptotic solutions are often not available. In those cases it is customary to transform 
the given boundary value problem, by a change of the independent variable, into a new 

boundary value problem which does not tend to a problem of lower order when p tends to 

infinity. 

as an example for such a transformation we take the differential equation (101), for 
the special case 

bg(x) > 0 . 

The case b Q (x) < 0 can be treated analogously. Without loss of generality we may further 
assume that 

a - 0 . (209) 

We shall refer to this boundary value problem as the problem (L). 

We now introduce the new independent variable 

s * (Sc , (210) 

and transform the boundary value problem (L) into an equivalent problem in x, to which we 
shall refer as the problem <L). Let U(x,p) be the solution of the problem (L). 0(x,p) 
or some of its derivatives will have a boundary layer at x « a » 0. The function U(x,p) 
is changed, by the transformation (210) into 

U(x, p) - U(£, p) . (211) 

P 

U(z,p) is the solution of the problem (iT). Since 

U <i) (x,p) - p _i U (1) (x,p) , (212) 
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the problem (L) can be written 


fir) 5 y (B> ♦ b 0 (2)y ( *) ♦ l e v (i)p-V n - V> ♦ l **» 

p v*l p 




W P 


(2.13) 


(X.) -X. 

y 1 (pg) - t lP 


0 < 1 < r 


L 1 S 


< V 

y (0) - t t p 


(214) 

r+1 < i < n . 

If we let p tend to Infinity In the coefficients of (213) we obtain the simple 
"limiting" differential equation 

y (n> + b„<0) y (B) - 0 . (215) 

It may be expected that the function U(x, p) tends with increasing p to a function 

u(s) - 11m <7(s,p) (216) 

p-M 


which satisfies the differential equation (215). Since (215) is of the same order as (213) 
we expect that the function V(x,p) will not have a boundary layer for large p. 

The transformation (210) may be described as a stretching of the function U(x,p). If 
U(s,p) does not have a boundary layer, we have, in a way, "stretched out" the boundary 
layer. 


The problem arises than what boundary conditions are satisfied by the limiting 
function u(z). If the boundary layer of the function U(x,p) at x ■ 0 occurs in 
U(x,p) itself and not in a derivative of U(x,p), the interpretation of (210) as a 
stretching which becomes infinite when p ♦ «•, suggests that we have 

u(«) - u<0) , 

where u(0) is the value assumed by the function 

u(x) - 11m U(x,p) (217) 

p*m 

at x « 0. u(0) will, in general, be different from the boundary value prescribed for 

U(x,p) at x • 0. 
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Tha transformation (210) is frequently used for tha solution of mors complicated 
boundary layer problems. Tha following points are than usually taken for granted without 
proof: 

(a) That lim U(z,p) exists. 

P-*** 

(b) That the limit u(z) satisfies the limiting differential equation 

(c) That u(«) - u(0). 

In our investigation we have been able to find an asymptotic approximation for 
0(x,p) directly, so that we did not need the transformation (210). But we are now able 
to prove the statements (a),(b), (c), for our problem (L). This is what we are going to do 
in this $. 


38. In chapter I we have derived for U(x,p) the following asymptotic representation 
(compare (171)): 

'V) P fo 

U(x.p) - (sn(X)Jp e ♦ (u(x)] . (218) 

In this formula we are using the following abbreviations: 


r w m \ \ Tr+ 1 *1 

* ’ f- b 0 < 0) > T 


(219) 


where 


i - (-D 


r+1 


vv • 

• • • • L,(u ) 

1 B 

VV • 

.... L r (U n ) 

WV 

• • • • wv 

w • 



( 220 ) 


(Compare formula (159)), 
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• • • » 


*1 !.,<«,) .... L,^) 

*2 W * * * * W 


t I. (U,) . . . . L (tt ) 
n n < n ■ 


( 221 ) 


(compare section 18); 


and 


*(x) - -b Q {x) 


( 222 ) 


u(x) 


m 5, 

V _!±1! 
L . 6 
irl 


« p (x) 


(223) 


(Compare section 17 for the 6^^, - we have dropped the Index b, used there, as 

unnecessary here -/ and compare section 6, theorem 1 for the u^{x)). 

In consequence of theorem 1, formulas (118) and (120), we may differentiate (218) 
formally at least n-t times, l.a. 


0 (1) (x,p) 


1-T 


• [* (x)n(x))p 


r+i 


£ *< 0*5 


[u U) (x)J, 


(1-0,1,...,n-1) 


(224) 


39. Knowing 0(x,p) and Its derivatives ws can now easily calculate the function U(r,p) 
and its derivatives with respect to s. For, considering that 

a £ /p *(C)d5 /? *<-£)di; 

a 0 - e 0 (225) 

we find, upon substitution of (212) into (224), that 


U ,l) (s.p) 


tiW> l (f) «n(*;)I P 
P P 


r+1 


f! 


p~ l < £ )], (1 - 0,1,... ,n-1) 
P 


( 226 ) 


Thus we have solved the problem (L). 
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The result of the passage to tha Halt In (226) can bs moat easily axprsssad by ons 


formula comprising all cases that can arise. If we Introduce the following symboli 

( 0. If t ^ 0 

e(t) - { (227) 

l 1, if t - 0 . 

Then we can write 

u(s) - 11a 0(a, p) - (i*W(0))e 2 ^ <0) e(T r( , 1 ) ♦ u(0) . (228) 

P*» 

Formula (228) is valid for 0 < s < «. 

We can now confirm the three unproved statements (a), (b), (c) of section 37. u(s) 
does exist, and it satisfies the limiting differential equation (215), as we may readily 
verify by substitution. In fact, u(s) reduces to a constant, unless 1 * 0. Finally 
we see from (228) that u( -) - u(0). This proves the stateaents (a), (b), (c). 

For the limit of the i-th derivative of U(x,p) we find, from (226), 

limO <i, (*.p) - U* l (0))e*'' ,<0> c(T rM ) , (1 • 1,2,...,n-1) , (229) 

p+m 1 

which is also the i-th derivative of u(s). We conclude from (228) and (229) that 0(0,p) 

J 4 \ 

and 0 (0,p) remain finite as p ♦ in other words, 0(*,p) does not have a boundary 

layer at x » 0, the boundary layer has been 'stretched out*. 

40. We have seen that the limiting differential equation (215) is satisfied by the 
function u(x) obtained by passing to the limit in the solution U(z,p) of the problem 
(L). How can we find a complete set of n boundary conditions satisfied by u(s)7 

If we formally let p tend to infinity in the boundary conditions (214), we obtain 
boundary conditions 



(-> - e(X i ) 


(i - 1,2 


r) 


(230) 











y 


(0) - e(T t ) 


(1 - r+1, r+2,...,n) 


(231) 


(T i ) 


From the remark made at the end of the laat section it follows immediately that u(z) 
satisfies the boundary conditions (231). If » a, the last boundary conditions (230) 
can only be satisfied if, by coincidence, i» u(0). We therefore replace the last 
boundary condition of (230) by the condition 

u(») - u(0) , 

which we have proved to be satisfied. The other boundary conditions (230) are certainly 
satisfied, since all the derivatives of u(s) vanish at z • *. 

$3. The Won-Homogeneous Differential Equation — S(y) ♦ M(y) ** f(x) . 

41. Introductions It is an open question whether the Main Theorem remains valid in full 

generality for the non-homogeneous differential equation 

- H(y) + M(y) - f(x) . (232) 

P 

But we will be able to answer this question in the affirmative when 

n - m - 1 , (233) 

provided none of the boundary conditions not canceled in application of the rule of the 
Main Theorem involves an order of differentiation greater than m-1. The meaning of this 
latter condition is easily understandablei If one of the uncanceled boundary conditions is 
of the order n-1, then the boundary value problem formed by the limiting differential 
equation 

M(y) - f(x) 

and the uncanceled boundary conditions is of a type to which the usual method of solution 
by means of the Green's Function cannot be applied, since this method presupposes that the 
boundary conditions are of lower order of differentiation than the differential equation. 

In our treatment of the homogeneous differential equation the relative order of 
differentiation of the limiting differential equation and of the remaining boundary 
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conditions did not play any important role, and it is by no moans certain that the 
condition above is really necessary in the non-homogeneous case. But it simplifies our 
proof greatly. 

42. The problem defined by (232) and (233) and the boundary conditions (108) will be 
called the problem (N). If, instead of (108), we prescribe the corresponding homogeneous 
boundary conditions, we shall speak of the problem (N'). We assume that 0 < r < n, i.e. 
we consider only actual boundary conditions leaving aside the Initial value problem. (The 
initial value problem can be treated by the same method.) It is easy to extend the proof 
below to the more general boundary conditions (201). 

Let z( x, p) be the solution of (W), *(x,p) the solution of (N') and U(x,p> the 
solution of the homogeneous differential equation (101) satisfying the non-homogeneous 
boundary conditions (108). Then 

Z(x,p) - z(x,p) + U(x,p) . (234) 

Since the asymptotic behavior of U(x,p) for large has already been investigated in 
chapter I it is sufficient to discuss the problem (N'). 

43. The Green's Function: 

It is known that the function z(x,p) can be written in the form 

z(x,p) » f® G(x,t, p) p f(t) dt . (235) 

Cl 

The "Green's Function" G(x,t,p) can be constructed in the following manner: 

Let y,(x,p), y„(x,p) • • • y (x,p) be a fundamental system of solutions of 
* 2 n 

— N(y) + M(y) « 0, and set 
P 
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Mt,p) 


y (n-1) (t p ) y (n . , . y ^ n ^(t.p) 

y‘ n_ 2 ) (t,p) y< n ' 2 ) (t,p) • • .y^ n_ 2 ) (t,p) 


y^t.p) y 2 (t,p) 


y n (t .p) 


k(x,t,p) - agn(x-t) 


y^x.p) 


y 2 (*»p> 


• • • y_<x»p> 

n 


yi n 2) (t,p) yi n 2) (t,p) • • • y* n 2) <t,p) 


y^t/p) y 2 (t,p) • • • y n (t,p) 


(236) 


and 


g(x,t,p) 


1 k(x,t,p) 

2 h(t,p) 


Than 


G(x,t,p) 


(- 1 ) 


H(x,t,p) 
A< p) 


(237) 


(238) 


whare 


and 


H(x,t,p) 


y^XfP) y 2 (x»P> 

W W 


y (x r p) g(x,t,p) 
n 

W (, (g) 


(239) 


w 


L (y,) 
n 2 


\ ( y„> 

n n 


L (g> 
n 
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J 


V y 1> V y 2 J • • * W 


l 2 (y,> i. 2 (y 2 ) • ■ • L 2 (y n ) 


l <y,) i* (y ) • • • l <y ) 
n i n 2 n n 


i “ 1 , 2 ,...,n, means that the operator 4 ia applied to g(x,t,p) considered a 


function of x. 


44. The asymptotic value of the Green's Function: 


We take as the fundamental system y^fx,p) the functions (114), (115) specialized for 


n-m ■ 1, i. e. 


p r v>(odc 

y^x.p) - v(x,p) - e (n(x)] 


y 1 + u <x,P ) " V* ,p) ” lu y <*)! 


(p - 1.2 .n- 1 ), 


*<x) - - b 0 (x) . 

(Note that our notations differ somewhat from those of chapter I and of chapter II, §1.) 


We have then 


L[(V), 


1 L i (k> 

V 9) " 2 h(t,p) 


4(0,). 


L l<Vi> 


(n- 2 ) ( n- 2 ) „(n _ 2 ) (a . . 

L (k) - V (t<p) °1 (t ' p) * • * Vi {t ' p) 


V(t,p) 


o,(t,p) ... a(t,p) 

I n 


If '^ f ’ >• > . 



4 
















with 


Tjet us assume that 


L i - L 1 • 

for 

i < r 

L*--Li , 

for 

1 > r 

h 0 (x) 

< 0 

* 


Then 


L’(V) 


X. X. 

P tv> 1 ( B) ti<6))e pw , for i < r 


T i T i 

. - p If (a) n< a) ] , for 1 > r 


where 


w - r B *<£)<*£ 

Cl 


If we expand (242) in terms of its first column, we see that in a < t < 8 the 
etc. terms of the expansion are of lower order than the second term. Hence 

1^00 - t-'(V)[h(t)l - V (n_2) (t,p) tl^OO) , 


where 


h(t) 


u< n - 2 »(t) u<"' 2> (t) • • • u (n ' 2) (t) 

I i n~* l 

u< n " 3) (t) u‘ n ' 3) (t) • • • u (n ‘ 3) (t) 

I 2 n -1 


U ? (t) 


u 2 (t) 


• • U (t> 

n-1 


and 


k(x,t) - sgn(x-t) 


u t (x) 


ui n - 3, (t) 


ll * • I U ( V ) 

2 n-1 1 ' 


(n-3), n-3 

J, (t> • • • u , (t) 

2 n-1 


u^t) 


u 2 (t) 


• * ‘ Vl (t) 


(243) 

3rd, 4th, 

(244) 

(245) 

(246) 
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Expanding h(t,p) in tarna of its first column we obtain similarly 

h(t,p) - V tn-1) (t,p) (h(t)], for a < t < 0 . (2 

h(t) i 1 0, because otherwise ^(t), Uj(t),...,u n _j(t) would be linearly dependent. 
From (237), (243), (244) and (247) we find 


L t (g(x,t,p)) 


2V ln-1) (t,p) 


g(x,t) 


X(x,t) 


2 h(t> 

Furthermore, expansion of (236) in terms of its first column yields 

k(x,t,p) - sgn(x-t) V(x,p) [h(t)) - V <n-2) (t,p) [k(t>] 
and therefore, because of (237) and (247) 


g(x,t.p) - sgn(x-t) —- —777 

2V (n1) (t,p) p {t) 


(g(Xft)i . 


Now we substitute (250) and (248) into (239) and see that 

H(x,t,p) - H 1 (x,t,p) + H 2 (x,t,p) 


H (x,t,p) - --- 

’ 2V (n - 1, (t,p) 


V(x,p) 

(u^(x)] » • • tu |l _ 1 (x)l, sgn(x-t) 

V(x,p) 

X.,(V) 

• • • 

(L 1 (u 1 >) • • • [L t (U n _ t )] 

x.;(v) 

• • 

• • • 

I. (V) 
n 

[t (U )I • * ■ [L (U )] 

n i n n - » 

• • 

I»' (V) 
n 
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V(X,P) 

lu^x)] 

• • • lu , (X) ] 
ft— 1 

Ig(x,t)) 


I» 1 (V) 

• • • 

Cl* 1 <u 1 ) ] 

• • • [L (u )1 

1 n— i 

tx. 1 (g)l 

• • • 

(253) 

• • • 

t (V) 

ft 

• • • • • 

(VVi 

l Wi» 

s s • 

(L n (i)] 



The asymptotic expression for A(p) Is (compare chapter X) 

X. 

A(p) - L 1 (V>[A] - (*,($> n(S) AJe pw 


W • * * ^2 (u n-1^ 


L (u ) • • • L (u ,) 
n I n n-l 


As In chapter I we make the assumption 

A * 0 . (256) 

(A is identical with ^ of formula (159) for this particular case.) 

45. Corresponding to the representation of H(x,t,p) as a sum of two terms in (251) we 
find, upon substitution of (251) into (238), 

G(x,tp) - G^(x,t,p) + G 2 (x,t,p) (257) 

with 


Gj(x,t,p) - (-1) 


n Hj(x,tp) 


H (x,t,p) 

n i 


G 2 (x,t,p) - (-1) 


Using (235) this leads to a representation of the solution s(x,p) as a sum of two 
integrals: 


-- ■* ■' •''-‘i. * .» 
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(260) 


with 


*(x,p) - z ^ (x, p) «• z 2 <x,p) 


z^x.p) - /® G 1 (x,tp) p f (t) dt (259) 

z,(x,p) - J B G (x,t, p) p f(t)dt . (260) 

4 Cl 2 

We shall prove that, as p tends to infinity, z^(x,p) tends to zero, while z 2 (x,p) 
tends to a solution of M(y) - f(x). 


46. In this section we are going to show that 

lira z. (x, p) - lira f B G,(x,t,p) p f(t)dt - 0 . (261) 

p-M» p-M» 


To this end we write (259) in the form 


z 1 (x,p) » /* G 1 (x,t,p) p f(t)dt ♦ / B G^(x,t,p) p f(t)dt 


(262) 


and prove that each of the two integrals in (262) tends to zero. 

a) In the first integral of (262) we have t < x. Hence, G^(x,t,p) has to be 
determined with + V(x,p) as the last terra of the first row of the determinant in (252). 

In order to find an asymptotic expression for H^(x,t,p) we expand the determinant in 
(252) in terras of the minors formed by its first and last columns. These minors are either 
zero or of the form 


± 2 iyv) L g (V> 


i 2 V(x, p)I« g (V) 


T < r 
e > r 


i.e. they are of the form 


\ +T 

. T e 


tql 


T 

e 

P 


Ui 

o |* f(5)dC 


fa' I 


(263) 


(264) 
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are 


where q and q* are certain numbers different from zero. Since and t^ + ^ 

greater than all the other X|s and t|s, it follows from (264) that, in a < x < B, the 
expression of highest order of magnitude among all the minor (263) is ± 2 L l( V) L^V). 
and no other minor has the same order. Hence, (252) can be written, - in a < x < t -, as 
follows: 


(x,t,p) 


T^-n+l 

(-l)"-^^ 1 I, .(VI 


-p f *{5>dc 


<f g4 ’ , (a) n(a) 

^"“’(t) n(t) 


v(x) 


(265 


where 


(u,(x)l.(u ,(x)] 

I ll" 1 

‘VV. l Wl J 


v(x) 


(L^tu, )1 

‘WV 

• • • 

(L (u.)J 

II 1 


Wi 11 


[t r + 2 < Vi n 


(L 


n (u n-i n 


(266) 


is the cofactor of the minor ±2 L^V) I^^tv) in (252). (265), (258) and (254) show 

that for t < x. 


G(x,t,p) 


(- 1 ) 


vr n+1 


-p /* <^(oae 


± ^(a) n(a) v(x) 

v> n_, (t) n(t) a 


(267) 


and therefore 


/* G(x,t,p)p f(t)dt 


T_.-n+2 

(- 1 ) r+1 * r+1 


(q) n(o) ^ ~ f *(e)dC P(t,p) 
A 


(268) 


with 

















F(t,p) 


tit) 

[* n_1 (t) n(t)J 


The integral in the right hand member of (268) tends to zero as p since its 

integrand tends to zero in the interior of the interval and remains bounded at the 

T r+l " n+2 

endpoints• 6 remains bounded, because we have assumed that 

T r +1 <n ' 2 ‘ 

Consequently, the left member of (268) tends to zero, as p ♦ 

b) In the second integral of (262) we have t > x. Hence, in the determinant in 
(252) the last term of the first row is -V(x, p). Expanding this determinant in terms of 
the minors formed by its first and last columns we see that this time the minor of highest 
order is formed by the two first rows of the determinant. The value of this minor is 

±2 V(x, p) L^v) . 

A calculation analogous to that used in part a) of this section leads to 

J B G (x,t, p) p f(t)dt - p ""* 2 [ A ' ** 2 ] (269) 

1 

where the constant Q is the limit of the cofactor of the minor above. Since t > x in 

the integral in the right member of (269), the integrand of that integral tends to zero in 

the interior of the interval of integration, x < t < B, as p approaches infinity. 
Furthermore, the integrand is bounded at the endpoints of the interval of integration. 
Consequently, the integral tends to zero. On the other hand, our assumption m > 0 

implies n-2 > 0, and the power of p in (269) is therefore not positive. Hence, the 

left member of (269) tends to zero, as p ♦ «. 

This completes the proof of (261), 
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47. Ne now turn to the asymptotic calculation of i 2 (x,p), (aaa (260)). No nota firat, 
on expanding tha determinant in (253) with reapect to ita firat column, that H 2 (x,t,p) 
can be written, in a < x < 0, in the form 


where 


1,(7) 


V x ' t ' p) * 7*0 [H<x * t)1 


(270) 


H(x,t) 


u, ( X ) 


W 


u n _ 1 (x) g(x,t) 


l 2 ( Vl> V«> 


VV 


L n (« ,) L(g) 
n n-» n 


(271) 


ie the limit of the cofactor of the element L.j(V) in (253). Substituting (270) in (259) 
it follows that 


Gj(x,t, p) 


(_t ) n i **<*>0 
P <^(t> *A 


and this formula, together with (260) givea 

v*,p) * <-u n f e ti^if(t)dt . 

*t)7i 


In this expression we may pass to the limit under the integral sign, since the asymptotic 
expression for the integrand is valid in the whole interval a < t < 0. Thus we obtain, 
replacing at the same time *t) by its value -bg(0. 


11m Z 2 <x,p) 

p*w 



H(x,t) f(t) 
A 


b Q (t) 


dt 


(272) 
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Now we combime (272) with (261) and (260), end conclude that the right member of (272) Is 
the value of lim z(x,t,p). But the definition of H(x,t) and A in (271) and (255), 

p-MO 

respectively, show us that 


n-1 H(x,t) 


(- 1 ) 


is exactly the Green's Function belonging to the differential expression 

1 


b Q (x) 


K(y) 


and to the boundary conditions ^ (i ” 2,3,...,n). 

Since these boundary conditions do not involve derivatives of higher than (n-2)nd 
order, this proves that z(x,p) tends to a solution of the limiting differential equation, 
satisfying all boundary conditions, except the first one. 


3 
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48. We summarize the results of thia ( In the following theorem 


Theorem 5i 

We consider the problem (N) defined by the differential equation 

- N(y) + M(y) - f(x) 232) 

P 

and the boundary conditions (108). N(y) and M(y) are differential 
expressions of the form (103) and (104) with 

n-m-1 . 

We make the following assumptions! 

(a) Conditions 1* - 5*, 7* and 8* of the Main Theorem (section 4) 
are satisfied. 

(b) b„(x) <0 in a < x < 8. 

<e> t r+f < n-1. 

(d) 0 < r < n. 

(e) f(x) is lntegrable in a < x < 0. 

Then the solution z(x,p) of the problem (N) tends - as p approaches 
infinity - to a solution of the differential equation 

M(y) - f(x) 

satisfying all the boundary conditions (108) except, in general, the first 
one, L,,<y). 

A strictly analogous theorem holds for (x) > 0. 
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Chapter III 


SOME RELATED PROBLEMS 

§1. An Example for Boundary Laver Problems In Systems of 
Differential Equating . 

49. A great number of unsolved boundary layer problems with Important applications can be 
formulated for systems of ordinary differential equations. He are going to discuss In this 
section a very elementary example In order to give an Idea of the boundary layer phenomena 
that can arise for systems! 

He shall discuss the system 


-1 „ 

p u" 

v" 


“ a u t b v 
• c u + d v 


(301) 


with constant a, b, c, A, assuming that 

a !< 0 . (302) 

As boundary conditions we prescribe 

u(a) ■ u , v( a) - v , u( B) ■ u_, v(B) - v„ (303) 

a a p p 

where u , v , u„, v„ are constants, 
a a 8 8 

The "limiting problem", obtained by setting p 1 - 0 in (301) is equivalent to the 
differential equation of second order 

v " “ a v - 0 (304) 


where 


and to the relation 




1305) 


( 306 ) 









The following questions ariset 

(a) Do the solutions U(x,p), V(x,p) of (301) and (303) converge as p + «? 

(b) What boundary conditions do the limit functions U(x) and V(x) satisfy if they 
exist, and are they solutions of the limiting differential equation (304)? 

The functions U(x) and V(x) cannot be expected to satisfy all four boundary 
conditions (303) and also the condition (306), for the prescribed boundary values (303) 
will, in general, not satisfy ths condition (306). The answer to these questions is 
supplied by the following Theorem! 


Theorem 6» 

Let u - U(x,p), v - V(x, p) be solutions of the system of 
differential equations 


p"V - a u ♦ b v >> 

v"»cu + dv-' 


(a, b, c, d constants) satisfying the boundary conditions 
u(a) » u # , v(a) - v a , u( B) “ u g» v < 


(301) 


(303) 


where u , u_, v , v_ are constants. Let us further assumes 
a B a 6 

Assumption 1: a f 0 . 

Assumption 2s a, b, c, d, a and B are given in such a way 
that the differential equation 


y" 



(307) 


with 

5 - I* 51 (305) 

and the boundary conditions 

y( u) » 0 , y< B) “ 0 

can be satisfied by the function y(x) = 0 only. 
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Then we can make the following statements: 

(A) The function V(x,p) converges, as p + ■», to a 
function V(x), satisfying the differential equation (307) 
and having the boundary values 

V( a) - v , V(8> - v, (308) 

a 8 

(B) If a > 0, then the function U(x,p) converges in 
a < x < 8, as p + to the function 

0(x) - - | V(x) (309) 


(which, of course, satisfies the differential equation (307)). 

(C) If a < 0, then the function U(x,p) does not converge 
but remains bounded, as p ♦ ", except when the prescribed boundary 
values satisfy the condition 


a u + b v « 0 
a a 

(310) 

« u B+ b v 8 - 0 , 


In which case statement (B) remains true for a < 0 also. 


50. To prove this theorem we start from the observation that U(x,p) and V(x,p) are 
both solutions of the differential equation 

p -, {y <4 * - <J y"} -ay'tiy-O . (311) 

In fact, if we multiply the second differential equation of (301) by b and substitute 
into it the expressions 

^ y p P ^ U . — ^ | J ^ y, ^ ^ ^ p ( ^ ^ ^ p W 

obtained from the first differential equation, we find the following differential equation 


in u alone 









I 


u -<pa + d)u"tpSu-0 , 

which is equivalent with (311). Similarly, it can be shown that (311) is satisfied by 
V(x,p), by eliminating u from the differential equations (301). It is also easily seen 
that the boundary conditions satisfied by U(x,p) and V(x,p), considered as solutions of 
(311), are, respectively. 


for 0(x,p) 


for V(x,p) 


y(a) « u 


y(B) - u„ 


y“(a) • (a u ♦ b v ) p 
a a 


y"(0) - (a u fl ♦ b v fl )p 


(312a) 


y(a) -v 


y(B) “ v. 


y" (a) - c u + d v 
a a 


y-($) - c U(J ♦ d v g . 


(312b) 


The boundary conditions for V(x,p) are of the type considered in our Main Theorem 

(chapter I, section 3). Applying the Main Theorem for these boundry conditions and for 

n-m - 2, bg(x) - -a, we see that for a > 0 as well as for a < 0 the function V(x,p) 

converges to that solution of the limiting differential equation (307) which satisfies the 

boundary conditions y(a) “ v , y(6) - v . This completes the proof for statement (A). 

a B 


51. For U(x,p) a special calculation is necessary, since the boundary conditions for 
U(x,p) contain p, a case not considered in our Main Theorem. Our principal tool in 
chapter I, the asymptotic representation of a fundamental system of solutions of the given 
differential equation, can be applied to the differential equation (311) and yields then 


(see theorem 1, section 6) 


U 1 (x,p) - e 


oVa(x-a) 


-o^a( x-a) 


Uj(x,p) - e 


Uj(x.p) - (A(x)l 


a <x,p) • (B< x)) 









Here the functions A(x), B(x) form a fundamental system of solutions of the limiting 
differential equation (307) and 


a m /p • 

With these functions we repeat the general reasoning of chapter X for this special case. 

We first re-write the boundary conditions (312a) in a form somewhat more convenient for our 
purpose: 


-1 

P 


y"( 6) 


1. “au„+bv„ 
1 8 8 


-1 

P 


y"( a) 


- a u + b v 
3 a a 


y( 8) - t 2 - u g 


y( a) * t. » u 

4 a 


Let us further introduce the abbreviations 

e^( 0-a) , ^ 


. e'^ ( a ' a) - W 2 . 


(314) 


(315) 


Plrst case : a > 0. 

Defining A( p) as in (134) we find 



(a)W t 

faJWj 

f'tA’lB)) 



A( p) - 

(1)W 1 

(1!W 2 

[A( 8)J 

[8( 8)) 

(316) 


[a] 

(a) 

p _1 tA-(o)] 

p'’[B"(a)] 



(11 

m 

[A ( a) 1 

[B ( a) ] 

• 

We expand this determinant in 

terms of 

its last two rows: 



A(p> - -[a 2 ] 

At 8) 

Ato) 

B( B) 

w » 

1 

B( a) 

-[a 2 ]D W t 

(317) 


with 


D 


A(6) B(8> 

A( a) B( a) 


In consequence of assumption 2 of the theorem to be proved we have 


( 318 ) 
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I 


Similarly, vie find 


4 1 

[a]W 2 

p _ 1 [v<e>i 

p" 1 [B*(0)l 

*2 

(i]w 2 

[a< en 

[B( 6)1 

*3 

[a] 

p _1 [A-(a)l 

p' 1 [B*(a)l 

*4 

HI 

[A( ail 

[B( a) 1 

taiw. 

S 

p" 1 [As (0)1 

P _1 [B«(0)1 

inn, 

*2 

[A( 0)1 

[B( 0)] 

[a] 

*3 

p 1 [A*(a) 1 

P _1 [B*(a)) 

[11 

*4 

(A< a)) 

[B( a)l 


-ta 1 1 D] 


- -[a * 3 D] 


In the determinant 


[alW, 

[aJW 2 

*1 

p" 1 [B"(0)l 

[ilw 1 

li]i» 2 

*2 

[B( 0)] 

[a] 

[a] 

*3 

p” 1 [B*( a) ] 

[11 

Ml 

*4 

[B( a) 1 


we subtract the a” fold of the first roe from the second row and the a fold of the 
third row from the fourth, thus obtaining 


[alW 1 

[a)W 2 

*1 

P _1 [B*(0)1 

[01W 1 

[0]w 2 

b 

- — V 

a a 

[B( 0) 1 

[a] 

[al 

S 

p _1 [B"(a)l 

[0] 

[0] 

b 

- — V 

a a 

[3(0)1 






Here we have used the fact that, in consequence of 013) 



Expanding (321) we find 


Similarly, we prove that 


Aj( 0) 


b 

a 



B< S) 

B( a) 


(a 2 )W, 


A<8> 

A( a) 



[a 2 ]W. 


(322) 


(323) 


(324) 


52. From (313), (317), (320) and (321) it follows that 


lim 

p-M» 


A,(p) 

A(p) 


O^X.p) 


0 


lim 

p-*°o 


* 2 (P) 

A(p) 


« 2 <x,p) 


0 


in o < x < 0 . 


(323) and (324) show then that 

U(x) - lim 0(x,p) 

P*» 

Is that solution of (307) which assumes the boundary values 

U(c) - - ~ v , 5(8) - 7 v . 

a a a b 

This proves statement (B) of theorem 6. 

Second case; a < 0 . 

If a < 0, then the absolute values of w, and W 2 oscillate with increasing p 
without tending to a limit. In this case we obtain from (316) the asympotlc expression 
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Co 


(325) 


A 2 ) (W, - V» 2 ) 

instead of <317). For A^p) and AjJp) we have 

A,(p) - -CO a] <Jt 1 - < 3 <* 2 ) 

Aj(p) - -CO a] (t 3 » 1 - t 3 ) 

and for A 3 <p) and A^(p) 


Aj (p) 


-C« 2 1 


b 

“ a V B 


b 

— v 
a a 


B( 6) 

B( a) 


(W, - w 2 ) 


a 4 (p) - -[a 1 


A( 0) 


A(a) 


b 

a V B 


b 

a v o 


<W. 


V 


53. From these expressions It follows immediately that 


4 A (p) 

U(x,p) - — 0 l (x,p) 


does not converge in this case. For 

77^7 <&3<P> ♦ & 4 (P) V*' p) > 

converges to the same solution of (307) as in the case a < 0, while the 


J, / • u um. i/T*l( x-a) VVS -io/TIl( 

A <4 1 U 1 + 4 2°2 a(W,^ 2 ) * + *^V^T) 


(326) 

(327) 

(328) 

(329) 


expression 

x-a) 
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does not converge unless 4^ “ t •• 0, in which case it is identically zero. 

This completes the proof of theorem 6. 

If assumption 2 of theorem 6 is not satisfied, then our reasoning does not hold any 
more. In that case it would be necessary to take into consideration also the second terms 
of the asymptotic series used, in order to find the order of magnitude of 4(p). 

Assumption 2 is easily seen to be equivalent, in this case, with 

Si>0 


and 


J- (f5-a) + Nrt when ~ < 0 


where N is any positive or negative integer. 

The case a » 0 could be easily treated by the same method. 

The boundary layers in this § occur only for the function U(x,p), which with 
increasing p tends to a function which does not have the prescribed boundary values, 
except, when these boundary values satisfy the condition (310). 

A more adequate and general treatment of boundary layer problems in systems of 
ordinary linear differential equations could probably be based on the asymptotic solution 
of linear systems as developed by Langer and G. D. Birkhoff [51. The assumptions of that 
theory would, however, have to be generalized for this purpose. 


$2. An Example for Boundary Layer Problems with 
Singularities in the Interior . 

54. Introduction! If the assumption 6* of the Main Theorem in chapter I is dropped, i.e. 
if we admit zeros of b Q (x) in a < x < B, then our whole theory becomes Invalid. For 
the zeros of bg(x) are usually singularities of the limiting differential equation 
M(y) - 0, and the theory of the asymptotic solution of differential equations, which was 
our main tool, fails in this case. 
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The general treatment of boundary layer problems In thla caee would probably require 
an entirely new approach. But It la already Interesting to Investigate a very simple 
special problem of this type in which the differential equation can be solved explicitly. 
We shall discuss the boundary layer problem of the differential equation 

- y" ♦ b(x)y' - f(x) (330) 

P 

with the boundary conditions 

y(o) - i 2 , y(S) « i, . (331) 

We make the following assumptions: 

1. b(x) is regular analytic in a < x < 0 . 

2. f(x) is regular analytic in a < x < 6 . 

3. b(x) has a zero at the interior point x « r of the interval a < x < B. But 

b'(r) 0. and b(x) does not have any other roots in a < x < 8* (This number r has, 

of course, nothing to do with the number r used in the first two chapters.) 

Assumptions 1 and 2 are by no means essential. We introduce them only in order to 
simplify our reasoning. 


55. The solution of the boundary value problem for the differential equation (330): To 
simplify the calculations it shall be assumed that the boundary values are 

y( a) - y( 6) - 0 . 

The case of non-hosK>geneous boundary conditions, which does not add any new features to the 
problem, la discussed in section 59. 

The general solution of (330), as obtained by elementary methods can be written in the 

form 


(J(x,p) - J* d? /* p f(n)e' p(M5) - Mn)) dn + c 1 J* e _pMe> de + c 2 


where 












A(x) - b(x) dx (333) 

and X, y, v, , cara arbitrary constants, which ara rastrlctad by tha prescribed 
boundary conditions. Tha five constants ara, of coursa, not assantial. In reality (332) 
depends only on two essential parameters, so that thraa of the five constants can be chosen 
arbitrarily. In order to obtain a form of tha solution suitable for the boundary layer 
problem. It Is convenient to set u “ 1, while the choice of X and v shall be left 
undecided for the moment. For typographical reasons It Is useful to Introduce the 
abbreviations 


p s “ £ de ^S 0 f <n)e" pU<C) " x(n)) dn (334) 

Q* - £ e‘ pX(C) dC . (335) 

Then (332) can be written 

0<x,p) - P* ♦ CjQ* ♦ c 2 , (336) 

and substitution of tha values x ” a and x » 8 into (336) leads to the two linear 
algebraic equations 


0 * p ® + c 0 ® c 
X 1 X 2 

0 ■ + e + e 

X 1 W X 2 


for c 1 and c 2 . Calculating c 2 from these equations and using the fact that 

Q® - q“ » Q®, one finds that 
A A a 


p x B ®x - p ?> 


(337) 
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investigation can ba obtained immediately fro* (337) by tha following conaidarationai 
Substitution of X for x in (336) shows that c 2 - u(X). As X was arbitrary, this ia 
true for any value of X, so that (337) can ba regarded aa the desired solution of (330) 
with X instead of x as independent variable. Writing x for X, the solution of the 
boundary value probla* is therefore obtained in the fora 


, °*„a «*„. 

a a 


56. The asyaptotic value of / F( x)e dx for large pi 

The solution of (336) is composed of integrals of the for* 

P(x)e pl * x> dx . (339 

It la therefore laportant to have asymptotic expressions for such integrals for large 
values of p. The following theorem, a proof of which can e.g. be found in a paper by 
O. Perron [8], will be the chief tool of the subsequent investigations. 

Theory*i If F(x) and <4x) are regular analytic in s < x < t, if 


C « 0, for x • R, a < R < t 
*<x) < 

v < 0, for x R, s < x < t 


<p(x) » (x-R) p (gQ + g^(x-R) + •••) 


is the Taylor series of (x) around x«R, then 


ft F(x)e pl * x) dx - if r(R)r4> - W) p _Vp , 

J S P P . |VP 


P l9 0 l 


where the brackets "( ]* have the meaning defined in section 5. (The number p here has 
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where tha brackets "t 1" hava tha meaning defined in section 5. (The number p hare hae 
nothing to do with the number p of chapter I.) If R " s or R-t, tha same is true, 
but with the factor 2 in (340) missing. 

This theorem can easily be generalized so as to include also the case ^(x) > 0: 
Theorem i If P(x) and *(x) are regular analytic in s < x < t, if *(R) is the 
maximum of *>(x) in s < x < t, where s < R < t, and if 

*>(x) - *<R) + (x-R) p (g Q + (x-R)g 1 * ...) 

is the Taylor series of ^(x) around x “ R, then 

f* F(x).^ lx) dx - [i P(R)r<i) —L—1 p- 1/P .*<*> . (341) 

p p lg 0 l 1/p 

If R » s or R - t, the same is true with the factor 2 in (341) missing. 

Proof i The integral 

£ P(x)e P< * <X) -* (R)) dx 

satisfies the assumptions of Perron's theorem with (x) - v (R) instead of ip(x). Since 

*(x) - *(R) - (x-R) p (g 0 + (x-R)g 1 ♦ ...) , 

(341) follows immediately, if (340), applied for the exponent <Mx) - f (R) is multiplied 
on both sides by e p . 

57. Passage to the limit in (338), if 


< 0 

for 

a 

< x < r 

b(x) - 0 

for 

X 

- r 

> 0 

for 

r 

< x < $ • 
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The asymptotic calculation of the integrals in (336) is easiest when the constant v 
in the definition of A(x), formula (333), is chosen equal to r. If thia is done, the 
function A(x) satisfies the following conditions > 


A(x) > 0, x j* r 
A(x) « 0, x - r 

A(x) is monotonic increasing for x > r 
A(x) is monotonic decreasing for x < r 


(342) 


M x) has therefore the shape indicated by the figure below. 



The passage to the limit, as p * «•, in (336) leads to different results, according as 
x is less, equal or greater than r. 

Case a). x > r . 

x B 

Application of formula (341) to Q and Q show that for x > r both integrals 

B a 

have the same asymptotic value, hence 


lim —^ » - lim —2 - -1 • (343) 

P+* 8„ P**° 

a a 


a 

To the inner Integral of formula (340) can be applied. For, in this case, 

r < x < n < 5 and, in this range, — <A<O - A( n)) as function of n assumes its maximum 
value 0 for n • £, on account of (342). As 


| n (-A( 5) + A( n) 


-b( n) 


is, by assumption, not zero, the number p of (340) is here equal to 1, hence 
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J* p f(n>e“ ptMO " M ’' M dii- (f(£) 


and therefore 


lim P® 

P-Mi X 


f® iiii 

j x b<£) 


*5 


(344) 


The expressions in brackets in (340) and (341) can be sure to be different from zero, 
if F(x) r 1 0 in a < x < 6* In order to avoid too lengthy formulas, this additional 
assumption shall temporarily be made. It is, however, by no means essential, and it will 
be shown later how to proceed without it. 

The letter E in the subsequent formulas shall be used to indicate non-vanishing 
constants. Note that the same letter E will be used for different constants. 

One finds immediately, by application of (341), 


Q 


B 

x 



-pMx) 


(345) 


in 

Q 8 - /».-** *>d C 
a 1 a 

the exponent reaches its maximum 0 for £ “ r. As, by assumption, b(r) * 0, but 
b'(r) f* 0, p is equal to 2 in this case and 

Q? - »I 7T • (346) 

a v p 

In order to find the order of magnitude of p“, consider that a < £ < n < x and r < x 
in the exponent -(A(£) - A(r\)) occurring in P°. Hence, the maximum of -(A( £) - A( n)) 
as function of n for fixed £ is 


-<A(£) 

- A( £)) - 0, 

for 

A(£> 

> A(x) 

”(A( £) 

- A(x)) > 0, 

for 

A( £) 

< A(x) 


The asymptotic value of the inner integral of P* is 
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/ e pf(n).- p(A( 5 , " MT ’ )) 4n 


IB] 


, for A( C) > A(x) 


[B}* -ptM5> " Ml0) , for A(£) < A(x) 


The contribution to of that part of the interval a < £ < x for which A( £) > A(x), 

(if it exists), can be neglected in comparison with the part where A(£) < A(x). As 


/“ .-P(A(5)-A(x))^ . m ^ 


1 pA(x) 


one haa therefore 


(B1 


1 pA( X ) 


75“ 


(347) 


From (345), (346) and (347) It follow*, finally, that 


P“ - IB] l 

Q 6 x P 

a 


and therefore, using (344), 


Um W(x,p) - - /® dC, for x > r 

P-H* 


(348) 


Caae b). x < r . 

A consideration analogous to that used for case a) leads to 


11 m U(x,p) - /* £ 7 ^ d£, for x < r . 

”(5) 


(349) 


See also case c) and d), section 58. 

Case c). x ■ r . 

Is regular analytic at x « r, then one shows Immediately, by a similar 
that 


if £rr 

b(x) 

consideration. 


2 b(f} ^ * 16 bff) d5 *' for x ” r 
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(The factor j is due to the Q's in (338).) in other words, U(x,p) tends in this case 
to the arithmetic mean of the two limits at x - r. 


58. Passage to the limit in (338), if 

C > 0, for a < x < r 

b(x) \ - 0, for x » r 

\ < 0, for r < x < 6 . 

In order to operate as much as possible with positive quantities it is convenient to 
set now 


A(x) - J* b(x) dx . 


(350) 


Then A(x) satisfies the conditions (342). If in the definition of P fc and o'" the sign 

8 s 

of the exponents is changed, U(x,p) can again be written in the form (338). 

In addition to the distinction between the cases x > r and x < r, the relative 

size of A( a) and A(8) plays now a part in the proof. 

Case a], x > r, A(8) > A(a) . 

I*t s > r be the value for x for which 


A(s) - A( a) 


Then 


(E] - e pA(X) , x > s 


[E] ’ e^ 3 ’ , x < 

0 


(351) 


and 


Q 8 - [E] - e pA<6) . 
a 0 

o 

The relative positions of x, £ end n in P® are indicated in the figure below. 


(352) 
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One sees that max(A(£) - A(r>)) - A<?) - A(x) for fixed i. Therefore 


/J p f(n>. pUU) - Mn)) dn- [■] e ptMe) - A(x)) 

rS b P<A<5)-A<x)) 2 e P (A( S)-A(X)) 

■*x P 


and 

_8 1 p( A( 8)-A(x)) 

P x - [E] - e 

Furthermore 

For p a consider again the relative position of x, £ and n* If t > 
Integral in P® remains finite, as p ♦ », hence only the case £ < r has 
considered# In that case 

max(A{£) - A(n)) - A(C) - Mr) - M C) 

and 



r p f(n) 


p(a< e)-R( h>) d 


tE] /p e 


pM£) 


As 


one has 


f° e pA<t) d t - 
J x 


r t E] -; e pA(x) 

\[I] - 1 p e pA(a) 


x > s 


x < 3 


(353) 

(354) 

the inner 
to be 
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V. 


1 pA( x) 

[El e , x > s 


[E] e , x < a 


Substituting formulas (351) - (338) ona obtains 


U(x,p) 


tE1 i 


[El A. a pA(X) 

*P 


(El - B p < A <a)-M*) > t _, 1 pA(a) 

1 J p " lE) /p 8 


lim U(x,p) - ±» . 

p-ooo 


Case b). x > r, A(B) < >(a) . 

The reasoning is of the same type as in case a), only the orders 
terms change. Let s < r be the value of x for which 


One finds 


*(a) - A( 8) 


—f, - [El . 

Q 8 

a 


Por p x the as y ra P totlc formula (353) holds unchanged. Furthermore 


Q S - [El i e pA(B) 

X p 


Q 8 - [El 1 e pA(a) 

a p 


The relative positions of x, £ and n in the inner integral of 

x 

figure below. Only the case £ < r has to be considered, as for £ : 


(355) 

x > s 

X < s 

(356) 

of magnitude of the 

(357) 

(358) 

(359) 

can be seen in the 
r the integral 


tends to zero. 
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For e < r 


jj P*(rt .<*M»-w»d,. m /p.^ n 


/ a . (A(5 ) d 5 - m ip<*<“> 
x p 




The expression* (353) and (357) - (360), Inserted In (338) lead to 


U(*.«) - - m Ae.^W 

O *0 


Therefore 


11a U(x,p) - *•• . 

P*~ 


Ca*e c) and d). x < r. 


One might repeat the preceding argument* in analogous fora and with the same final 
result for x < r. Instead, one can also proceed a* followsi The transformation 


t m r - X 


changes the differential equation (330) into 


+ * (i) £■ ?(z) 

02 


b(*) » -b(r-s) 


f(*) • f(r-s) 
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The boundary conditions yta) * y( B) “ 0 are transformed into 


As 


y(r- B) “ y( r-a) - 0 . 


f * 

0 , 

for 

r- B < * 

< 0 

b(z) \ » 

0 , 

for 

z 

- 0 

l < 

0 , 

for 

0 < z < 

r-a 


The results of the preceding section can be applied to (362) and lead to 

11 m U(z,p) • ±«* , for 0 < z < r-a • 

P-m* 

But this is equivalent with 

lim U(x, p) - ±“ , for a < x < r . 

P-h. 

This method might also have been used to prove formula (349). 


59. Some minor generalizations: 

1) If f(x) is allowed to have roots in a < x < B» the asymptotic values used for 

p and P p are not correct for these values of x which are roots of f(x). If, in 
x x 

particular, f(r) ” 0, some of the expressions would even be incorrect for all x. The 
reason is that f(x) or f(r) appear in the constant factors occurring in these 
expressions. But in that case the calculation could be carried through with a slightly 
more general form of formula (340) also contained in Perron's general formula in the paper 
( 8 ] quoted above. The result is again the same. 

2) Non-homogeneous boundary conditions: If to (338) is added a solution of the 
homogeneous differential equation 

- y" ♦ b(x)y' - 0 (363) 

P 

satisfying the non-homogeneous boundary conditions 

y(o) - tj, y(S> - l, (364) 

one obtains the solution of (330) satisfying the boundary conditions (364). 
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4 

Mr 

* 


As tha transformation 


y - y - t. 


changes (330) into a differential aquation of the same type and transforms the non- 
homogeneous left hand boundary condition Into a homogeneous one, It Is no loss of 


generality to assume that 


s-° 


The general solution of (363) can be written In the fora 


U(x,p) • c ( Q* + c 2 


The given boundary conditions lead then to 


°1 “ li ' c 2-° ' 

w a 


0(x,p) - — t . 

Q® 1 


To (365) the methods of sections 57 and 58 can be immediately applied with the following 


resulti 


r < 0, x < r 
l > 0, x > r , 


0 for x < r 


f 1. for x > r 
lim tJ(x,p) “ < 

p-M« v. 0 for x < r 

(> Q, for x < r 
b). b(x) ( 

l < 0, for x > r , 


a.) A( B) » A(a) , 


U(x,p) - 


[«] .P<* ( x>-A(B)\ x > , 
[«1 .p<*<«>-a(B)) # x < , 



.« • . 









M 0) < M a) 


6 .) 

lia 0(x,p) - t 1 . 

P*** 


60. Theorem 7 i 

Given the differential equation 

iy* ♦ b(x)y' - f(x) (330) 

P 

where b(x) and f(x) are regular analytic functione in the 
interval a < x < 0 and b(x) has exactly one root x - r 
in the interior of the interval, while b'(r) 4 0. Then the 
behavior for great valuea of p of the aolution of (330) which 
aatisfies the boundary conditiona 

y(a) - t 2 . y<0) - i, (331) 

depends essentially on the shape of b(x): 

C < 0 for a < * < r 
1). If b(x) < 

' > 0 for r < x < 0 , 

then the solution n(x,p) converges with Increasing p in 
the whole Interval a < x < 0, except possibly at x * r, the 
limiting function being composed of the two solutions of the 
limiting differential equation of the first order 

b(x)y' - f(x) 

satisfying one of the two prescribed boundary conditions. If 
these functions are bounded at x « r, then the solution of 
(330) converges at x » r to the arithmetic mean of the two 
limits at this point. 

C > 0 for a < x < r 
2). If b(x) < 

v < 0 for r < x < 0 , 
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* 


then D(*,o) dlvargaa with lncraaalng p at all pointa of tha 
lntarval. 
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Appendix 

A SHORT REPORT ON THE ASYMPTOTIC SOLUTION OP 
ORDINARY DIFFERENTIAL EQUATIONS INVOLVING A PARAMETER 

61. The aaln mathematical tool used In thia paper la the theory of the asymptotic solution 

of differential equatins Involving a parameter P for large values of this parameter. The 

most important results of this theory are contained in papers by G. □. Birkhoff [1], 

Noaillon [2], Tamarkin 13], [4], Langer [5], and Trjitsinsky [61. The asymptotic 

developments used by Birkhoff and Tamarkin. although of a very general character, do not 

apply to the particular differential equations of this investigation, because they assume 

that a certain "characteristic* algebraic equation formed with the coefficients of the 

differential equation has no multiple roots, an assumption not satisfied in our case. 

It would probably not be difficult to modify the methods used by Birkhoff and Tamarkin 

in such a way that they cover our case. But this is not necessary, since the type of 

differential equations considered by Noaillon and Trjitsinsky includes the differential 
equation (101). 

A complete proof of the main theorem of Noaillon and Trjitsinsky would be beyond the 
scope of this investigation, even if we restricted ourselves to the special case in which 
we are interested. We intend here only to give a summarised report on the methods of this 
theory and to show how, assuming the theorems proved, the asymptotic expressions of theorem 
1 can be obtained in our case. 

The theory consists of two parts. In the first part, which may be called the "formal 
part", the "exact" differential equation 

L(y,p) - 0 (401) 

which is essentially equivalent to the differential equation 

— ^l p) - - 0 (402) 

y 

is replaced by the "asymptotic" differential equation 

L(y,P) „ 

Y j 


0 


(403) 











Her* the symbol "j" has the following meaning] 



'1 

• * 

Jt 


f(x,p) J g(x,p) , 

where j Is a real number, stands for 

f(x,P) - g(x,P) - . 

p" 3 

As previously, E(x,P) Is a function such that there is a positive real number R so that 
|R(x,P)| is uniformly bounded for a < x < 0 and |p| > R. Usually, f(x,P) will be 
regarded as asymtotically equal to g(x,P) only if j is positive. But sometimes, e.g. 
in the case of the differential equation (101), a solution of (403) with a negative j 
will be an asymtotic approximation in the ordinary sense of the exact diferentlal equation 
(401). 

It is then shown that under very general assumptions a function Y(x,p) can be 
constructed which satisfies the condition (403), provided the number j is not too large. 

If certain differentiability conditions are satisfied, j may have an arbitrarily 
large value. This is the case treated in detail by Noalllon and Trjitxinsky. The case in 
which there is an upper limit for j is only mentioned occasionally by these authors. But 
since we are only interested in the first terms of the resulting asymptotic expansions, it 
is unnecessary to assume indefinite differentiability of the coefficients of the 
differential equation. This assumption is required only, if we are interested in the 
unlimited asymptotic expansion. Going over Noaillon's proof it is easily seen that 
assumption 4* of Noaillon's theorem in section 62 below is sufficient to guarantee the 
existence of the first term of the asymptotic solution. 

In the second, the "functional'' part it is shown that the solutions of the asymptotic 
differential equation (403) are asymptotically equal to the solutions of the exact 
differential equation (401). 

Essential for our application is furthermore the result that the derivatives of these 
asymptotic solutions of (401) are asymptotically equal to the derivatives of the 
corresponding exact solutions. 
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Th« complete etateaent of the result* of Bouillon's paper, as far as they are 


Important for our purpose, follows: 


62. Soalllon's Theorem . 

Part Ii Given Is the differential equation 

WfiP) “ l VAx, p)y <n " 1) - 0 (404) 

1-0 

satisfying the following conditions: 

1*. p is a real positive parameter. 

2*. x is a real variable 

3«. In the domain a < x < 8, p > R < a, 8, R constants) the coefficients P^fx,p) can 
be expanded in convergent series of the fora 

H — 

P A (x,p) - p 1 l B (x)p“* . (405) 

s-0 

(The are positive Integers.) 

4". The functions B 1# have at least n continuous derivatives in a < x < 8. 

5«, The coefficients B ao (x) in (405) does not vanish in any point of the interval 
a < x < 8* 

To these conditions 1* - 5* a sixth assumption has to be added, which can be most 
easily defined in the course of the construction of the asymptotic solution. 

If these conditions are satisfied, then there can be constructed solutions Y(x,p) of 
(403), each of which can be Witten in the form 

Y(x,p) - T*u (406) 

where the "principal factor" T is a function of the form 


i 

with 


T 





(407) 


- f A (X)p 01 
i-1 


i Hx,p) 


(408) 
















1 


w&mm 


the q. bain? non-negative decreasing rational numbers and p being a positive integer 
independent of j. u stands for tha "secondary factor" 

j. . 

u(x,p) - l * (x)p M • <409) 

v*0 

Here M is a positive Integer independent of j. The positive integer j' depends on 

j and increases with }. (In the application to differential equation (101) the first 

term of (409) is obtained if j » - -, as we shall see.) 

n-« 

As we have said before, the maximal value of j for which a solution Y(x,p) of 
(403) can be constructed depends on the number of times the coefficients B ig can be 
differentiated. It can be only determined in the course of the successive construction of 
the terms of (409). If the B ig can be differentiated indefinitely, then j and j' can 
have arbitrarily large values. 

Remarki The theory remains valid if the series (40S) are not convergent but only 
asymptotic expansions, nut we do not need this case for our application. 
functional Part : bet Y(x,p) » T«u be a solution of (403). Then there la a solution 
y(x,p) of the exact differential equation (401) such that 

- 11 

y(x,p) - T(u + p M [0]) (410) 

and this equation can be formally differentiated at least n-1 times, i.e. it can be 
differentiated treating the symbol [0] as if it were a constant. 


63. The construction of the asymptotic solution of the differential equation: 
A) The principal factor . 

The first step: In ^ substitute an expression of the form 


(compare (407), (408), (409)). The result of this substitution is an expression of the 


* T 1 \ 
l f (X)P ♦ P *[0] 
i-1 


4 

4 

*4 



i... v-. - vV ’$u. 


If 


* -v>. 










a 

the first term A^(x)P 1 in (408) in such a way that the coefficients fj(x) of the term 
of highest order in (411) vanishes* 

Noaillon gives a general method that allows us to determine the exponent and the 

function A^x) systematically in the general case. But our application being of such a 
simple type, and A^(x) can be found, in our case, more directly. 

Second step: Having chosen the first term A^P of ♦ we set 


P 


y 




where is the denominator of the rational number . This substitution transforms the 


expression 


into an expression 


w 


We now repeat the reasoning of the first step with 


L 2 (y 2 > 


instead of ** ** * >, in order to 

y 


find the term A P , considering, of course, for a only values that are smaller than 
2 2 a 

a. we continue in this manner until we arrive at the last term A P P for which a is 
1 P P 

still positive. That there is such a term, i.e. that we always attain an exponent 
a . < 0 after a finite number of steps, is proved in NoailIon's paper. 

p+l 

B) The secondary factor . 

First step: In L(Y) substitute Y - T'u where T is the function calculated in A) 

and u as yet undetermined. The result of the substitution is an expression of the form 

T K(u,P), where K(u,P) is a linear differential expression in u whose coefficients are 
1 /M 

power series in o « P , M being the common denominator of all the exponents in 

♦. In K(u,P) collect the terms of highest order in o. Then K(u,P) can be written in 


the form 
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K(u,P) - S G(u) + H(U/P) 


(412) 


Here S ia the highest power of o occurring in K(u,P) and G(u) and H(u,P) are 
differential expressions, G(u) being independent of P. It can be proved that the 
highest order of differentiation occurring in G(u) ia greater than tero. 

Second step: Find a solution of the differential equation 

G(u) - 0 . (413) 

We take this solution as the first term Eq(x) of the series (409). Since we want Bg(x) 
to be bounded in the whole interval of x in which we consider the asymptotic expansion, 
we have to add to the assumptions 1* - 5» in section 62 the condition 6*: The coefficient 
of the highest derivative in (413) does not vanish in any part of the interval ot < x < 6. 

Third step: In order to find the function E^(x), W > 0 of (409), Noaillon proceeds 
as follows: He determines by recursion a sequence of functions w^(x, ) wj(x, ),... from 
the formula 

G(w v > - | H(w v _i # P) . (414) 

It is easily seen that the w^tx.P) are of the form 

W - I s " 1 g v .(x) . (415) 

i-0 

The functions E^lx) are then given by 

P 

e u (x) “ I 9 VU (*> . (416) 

v* i 

It is not difficult to prove that the function (406), if determined by the 
construction which we have just outlined here, does in fact satisfy the relation (403). 

The construction of the asymptotic solution Y(x,P) is by no means uniquely 
determined, and it can be proved that the construction yields asymptotic expansions which 
are asymptotic approximations to a fundamental system of solutions of (401). 
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64. Asymptotic Solution of L(y) = N(y) + M(y) - 0. 

A) The principal factor. 

First step: We now apply the construction described in section 63 to the particular 
type of differential equation under consideration. Substituting 

y » Y » T»u 


MXi£i ( where L(y,p) is now the differential expression (101), and T and u are 

expressions of the form described by (407), (408), (409), we see that, unless i|i « 0, the 

term of highest order of ^- is (Aj(x)p ) . The condition that the highest terms of 

y N( y) M( y) 

-as function of p, and of —*— cancel out is therefore 

P Y Y 

a.n-1 a.m 

a"(x)p h -b Q (x) a“(x)p , 


o^n - 1 - a^ 


(417) 


and 


» n V. 

1 “ "Vi ' 


(418) 


provided iji / 0. From (417) and (418) we conclude 

1 


1 n-a 


(419) 


and 


A, - (-b 0 ) 1/n_n 


(420) 


Second step: Now we write 


o ft 


y =» e 


(421) 


where o is defined by p 3 a and (x) is one of the functions y (x) defined in 
theorem 1, section 6. Yj is the function 


ft MOdC 

0 

y 2 - e • u 


Here u is again the series (409) and 

P 

<M*> - l A (x) p 
i-2 


(422) 


(423) 
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He find 


v (v) „ o W - 1 ^ 1 (x)y- 

1 - * o (x) + V — + • • • • 

y y. 


(424) 


The term of highest order in the second term of the right member, if Y^ is replaced by 
its value according to (422) and (423) has the value 

, v-1 °2 v-1 

v A 2 a p if 

V V 

This term is of lower order in p than the term a <t , since 

*(»-»)«, 

p o ■ o -1 


and 


in consequence of the assumption 


(n-m) ol - 1 < 0 , 


^ < °t " n-m • 

The omitted tens in (424), which are indicated by dots, are similarly seen to be of still 

b(yj) Inn mm 

lower order. The terms of highest order in -, i.e. — a and b.o « cancel out, 

y 2 P 0 

since i(x) and o have been chosen such as to achieve just this. The next terms are 

(425) 


1 n-1 °2 n-1 

n p V p * . tor - y 


1 “‘V 


and 


„ „ -1 “2 -1 , M<y 2 > 

bA o p <e , for —- 

u * y*i 


(426) 


provided <jj > 0. (For, if “ 0, there are more terms of the same order as (425) and 
(426).) 

Following Noaillon's construction we try to determine and A 2 in such a way that 

these two terms cancel out. But setting the sum of (425) and (426) equal to zero, and 
inserting for *(x) its value (-bg(x) ) 1 ' /n “® leads to 

n-m - 0 , 

which was excluded. 

Hence, cij > 0 is impossible and therefore • <w(x). 
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B) The secondary tern. 


He follow the construction of Noaillon In order to find the differential expression 
K(u,p) of (412). Since 

(v) v v . vl v-1 , . 

y ■ Tlu «ii ♦ v« v u* •»■••• •) 

where the dots Indicate terns of lower order, we find by an easy calculation 


L(Tu) • T K(u,p) “ T{o m ^(a^ 11 ^u + nf” 'u* + b^® ^u + nb^ B ” 1 u') + H(u,p)} 


Hence, the function Bg(x) In this case la a solution of the differential equation 
obtained by setting the factor of o equal to zero. This differential equation can be 

written 

(a^ ♦ b^Ju * (np ♦ rab^u' 

or 

(a, b 0 - b,)u » b 0 (n - m)u’ . 

Therefore 

f x V g)b o <C) * V S) ir 

" *o b 0 <5)(n-n) ^ 

B 0 (x) - n(x) - e . (427) 

Interested In the other terras I^x), v > 0 of (409). 

1* - 6* of Noalllon's theory are satisfied In our application for the whole 
< 8. Condition 5», in particular, is equivalent to condition 6» of the 
chapter I. Hence, we conclude that there are n-m solutions of (101) of 

o F * <«« 

Y v (x,p) - e a v (n(x)l, *(v - 1,2,...,n-m) . (42S) 

Note that the function n(x) Is the same for all Y y (x, ) and that It does not vanish In 
a < x < 6. 


He are not 
Conditions 
interval a < x 
Main Theorem of 
the form 
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But w« can find more asymptotic solution* of th« dlffarsntlal aquation (101) by 
dropping th* assumption that th* i(i(x,p) of (408) is not ssro. In fact. If th# principal 
factor T of (40?) is equal to 1, the method used for the construction of th* secondary 
factor In section 63, B) leads to asymptotic solutions given by series of the fora (409). 
The first term of each of these series is a solution of the differential equation 
M(y) - 0. Taking a fundamental system u^(x) of m Independent solutions of this 
differential equation as first terms of m asymptotic solutions of (101) we can add to the 
n-m solutions of (101) given by (428) a more solutions of the form 

V(x,p) - lu (x)} , (v- 1,2,...,m) . (429) 

n-a+v 

linerly Independent, for sufficiently large p, can be given by calculating the asymptotic 
value of the Wronskian of these n functions. It does not offer any difficulties. 

This finishes the proof of theorem 1 of section 6. 
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